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Abstract 

In recent years, constructive field techniques and the method of renor- 
malization group around extended singularities have been applied to the 
weak coupling regime of the Anderson Model. It has allowed to clarify 
the relationship between this model and the theory of random matrices. 
We review this situation and the current program to analyze in detail the 
density of states and Green's functions of this model using the supersym- 
metric formalism. 



1 Introduction 

This small review is devoted to the elementary theory of random matrices and 
to the link between this theory and the Anderson model of localization/diffusion 
of a quantum particle in a random potential. 

More precisely we recall first the basic result of random matrix theory, 
namely the Wigner's semi-circle law for density of states, and give its rigor- 
ous supersymmetric derivation. 

Then we review the Anderson model, introducing the phase space approach 
to this model pioneered by Gilles Poirot, and summarizing the results of pQ. 

Finally in the last part we propose some generalizations of the flip random 
matrix model of [I] which are closer to the real Anderson model, using in par- 
ticular some hierarchical approximations. The control of such more realistic 
random matrices models is an important step towards rigorous theorems about 
the Anderson model in the weak potential phase. 

We thank Abdelmalek Abdcsselam, Jean Bellissard, Jacques Magnen and 
Gilles Poirot for collaboration and discussions on various aspects of this work. 
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V. Rivasseau also thanks F. Klopp for an invitation to lecture on this subject 
at Paris XIII- Villetanneuse University at the origin of this review. 



2 Random Matrices and Wigner's law 
2.1 The GUE 

The simplest ensemble of random matrices is the Gaussian unitary ensemble. 
It is a probability measure on random hermitian N x N complex matrices. 
Each coefficient in the upper triangle of the matrix is identically independently 
Gaussian distributed. Here the matrix is H = Hij, H = H* , hence Hij = Hji, 
and 

1 N 

P(H) = -exp(--TrH*H) (2.1) 

Z being a normalization factor. The matrix H is made therefore of N(N — l)/2 
complex variables Hj with i < j and N real ones Ha, so there are N 2 real 
random variables in H. 
Since 

N N 

TrH*H = ^ \ H V I' = E H l + 2 E( Rff « + ) ( 2 - 2 ) 

i—1 j — 1 i i<j 

we have 

Z = 2 n ' 2 (tt/N) n2 I 2 (2.3) 

and the covariancc rule is 

<HijH*, r >=^5u,5jj,. (2.4) 

The scaling factor has been chosen to keep the typical eigenvalues of H of 
size O(l) as N — > oo; indeed the typical size of the eigenvalues of a random 
matrix with covariance 1 is obviously of order \fN, by the law of large numbers. 

Physicists would like to know the statistics of the eigenvalues of H, and they 
are particularly interested in the two first moments of their distribution, called 
the density of states and the twodevel correlation function. 

The density of states, v{E) for a Hermitian matrix H is the quantity which, 
when integrated from — oo to A, counts the number of eigenvalues of H which 
are lower or equal to A. Since H has exactly N real eigenvalues Ai, Ajv, we 
have 

1/(25) = ^TrS(E H) (2.5) 

so that v{E)dE = 1. We can use the standard formula for the Dirac 

distribution 

S(x-a) = -- lim 3? . (2.6) 

7r £->o+ x - a + ie 
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Hence 

v{E) = — lim 3Tr- -. (2.7) 

v ' irN e^o + E - H + ie y ' 

Physicists call (E — H ± ie)^ 1 respectively the retarded and advanced Green's 
functions for the Hamiltonian H. 

The averaged density of states < v{E) > is therefore 

< v(E) >= - l_im / W^M^, (2.8) 

and < v(E) > dE clearly represents the probability for an eigenvalue of H to 
lie between E and E + dE, with normalization condition J < v(E) > dE = 1. 
The main results on the GUE ensemble is Wigner's semi-circle law: 

lim < v(E) >= ^^Vl-£ 2 /4 . (2.9) 

jV— «x> 7T 

The corresponding curve is really a semi-ellipse, but of course could be changed 
into a circle through a slight reparametrization of the covariance of H. The 
normalization taken here corresponds to J E 2 v(E)dE = 1. 

Wigner's law is a central result. It has been called the non-commutative 
analog of the Gaussian law of large numbers [2], and has been proved to hold in 
much more general cases than the GUE, for instance for band random matrices 

The next quantity of interest is the 2-level correlation, which allows to know 
the conditional probability to find an eigenvalue of H near E knowing already 
that one eigenvalue sits at E. More precisely it gives the probability to have two 
eigenvalues separated by an interval of width ui centered at E, and is therefore 

_ <-tf-y t .ffl> (2 , 0) 

In the GUE, eigenvalues are not independent but tend to "repel" each other. 
This is seen in the following behavior of the 2-level correlation i?2 

■ 2 

cin 7T H 

Ra(s)=5( 8 ) + 1- 1 ~, (2.11) 

where s = w/A, and A is the mean level spacing A = 1/iV < v{E) >. The 
delta function simply expresses the constraint of presence of an eigenvalue at 
E. Independence of the eigenvalues would mean lim s ^o R2(s) — S(s) — 1, hence 
no change in the probability for a second value to sit near E if a first is present. 
But here we have lim s ^o -Ra(s) ~ S( s ) — because of the s ^" s ^ s term. Hence 
there is chance for a second eigenvalue to sit near E if a first one sits at E. 
This is the phenomenon of " eigenvalue repulsion" . 

Physicists got intuition of this repulsion by the simple observation of the 
Vandermonde determinant that appears in the Jacobian of the transformation 
from the initial coefficients of the matrix to the diagonal eigenvalues and the 
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unitary diagonalizing matrix. In rough terms, we can diagonalize an Hermitian 
matrix H through a unitary matrix U : 

H = UAU* , U* = U~ x . (2.12) 

Then one can write the initial measure P{H)dH in terms of the coefficients 
of A and U . Clearly the measure on the unitary group must factorize from 
the eigenvalues measure since P{H) is invariant through action of the unitary 
group. Let us explain by a simple argument the well known result 

N 

P{H)dH = d/i(l/)e-T ££i a * Y[(\i - Aj) 2 JJ d\ . (2.13) 

i<j i—1 

To understand the appearance of the non-trivial Vandermonde factor rii<j(^ — 
\j) 2 (in addition to the ordinary trivial factor P(H) = e~Tr^i=i A *) we need 
only to compute the Jacobian at origin from the H variables to the A variables 
and the variables parameterizing U near the origin. For this purpose, we can 
derive the relation U* U = 1 with respect to a set of local parameters U r for a 
local chart of the unitary group near the origin. This gives 

S r = U*§- = -S* r , (2.14) 

the tangent space to the unitary group at the origin being the anti-hermitian 
matrices. From H = UAU* one finds 

ou r = ou r AU +UA ou;> (2 ' 15) 



hence 



Furthermore 



{U*—U) t3 = (S r A - AS r ) tf = (S r )ij(Xj - Ai). (2.16) 



so that the Jacobian to compute is 



J = 



dU r dU r 



1 

f(\,U) (\j — Xi)S r (U) 



g(U) n(A,-A,) 2 (2.18) 

i<j 



Clearly the presence of this Vandermonde determinant means that the eigen- 
values of a random matrix in the GUE case are not independent, but repel each 
other since the measure vanish at coinciding eigenvalues. Physically this level 
repulsion is analogous to some kind of Pauli exclusion principle between eigen- 
values, or to some two body logarithmic interaction: 

N N 
e -f Ef =1 ^ J| (A . _ A .)2 -Q ^ = e -f E f =1 A?+2 Ei< , loglA.-A,! J| ^ {2AQ) 

i<j i—1 i—1 
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which is analogous to Coulomb repulsion in two dimensions (also logarithmic) . 

It is possible to use the theory of orthogonal polynomials to analyze the 
large N limit and recover Wigner's law for this system or for more complicated 
non-Gaussian measures on H (for the GUE, orthogonal polynomials are simply 
Hermite polynomials). This is e.g. done in [J. See also [5] for another reference 
book on the subject. 

In this lecture we prefer to stress the supersymmetric approach to this prob- 
lem. It makes particularly transparent how Wigner's law results from a mean- 
field theory and a saddle point expansion which expresses the subtle balance 
between the Gaussian and Vandermonde terms in 12. 191 



2.2 Supermathematics 

In this presentation we follow the excellent concise review by Mirlin [Sj. 

Grassmann or anticommuting or Fermionic variables are pairs of independent 
variables, which for convenience are noted as complex conjugates, xii ■■■■>Xn> 
Xi, — ,Xn w ith the following properties 

XiXj = -XjXi , XiXj = -XjX* , XiXj = -XjX* , (2-20) 



dXi = J d X * = 0, j XidXi = J X*d X * = (2.21) 

Conjugation for Fermionic variables is not involutive but antiinvolutive, and it 
does not reverse the ordering of a product: : 

(X*)* = -X , W)*=XY ■ (2-22) 
Then for any N x N matrix M we find using these rules that 

J lidx*d Xi e-x* M x = det(^), (2.23) 

where x is now considered a vector {xti —>Xn)- For ordinary, also called bosonic, 
complex conjugate variables S, S* we would have got instead, requiring 5RM > 
for convergence: 

J XldStdSie- 3 '** 8 = det- 1 ^). (2.24) 

i 

Supersymmetric expressions are obtained by assembling bosonic and Fermionic 
variables in a symmetric way. For instance a supervector is 

$=(S u ....,S n ,xi,-Xn) , $* = (Sl....,S* N ,X*i,-X%)- (2.25) 
A supermatrix is M = I a , J , 27V by 2N, in which a and b are ordinary 



P b 

bosonic and a, p are anticommuting variables. We use Latin indices from 1 to N 
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and Greek indices with values b and / to distinguish the bosonic and Fermionic 
parts, so that we write 

M = M a /3 t ij , M b b,ij = a,ij , Mffjj — bij , M b fjj = ay , Mf^ij — Pi 3 ■ (2.26) 

Traces and determinants generalize into supertraces and superdeterminants: 

StrAf = Tra - Tib , SdetM = det(a - ab^p). det b~ l . (2.27) 
We still have 

Str log M = log SdetM . (2.28) 
Furthermore for any supermatrix (requiring tflMbb > for convergence): 

J d$*rf$e-** M * = SdetAf- 1 . (2.29) 

Remark that for a supersymmetric supermatrix M, i.e. one in which p = a = 
and a = b, this integral is normalized, namely SdetAf = 1. 
The inverse of a supermatrix can be computed as 

(o-<rfrV)- x -{a-ab- l p)- l ab- 1 \ 

-&- x p(a - ab- l p)~ l b-\l+ p{a-ob- l p)- l ob- 1 ) ) ' l ^ UJ 

Like for an ordinary matrix we can define a resolvent, or two point function 
J *ai*« e -** M * = (A/" 1 )?/ SdetM" 1 . (2.31) 

Further references on supercalculus, supermanifolds can be found in [7j. 




2.3 Wigner's law 

Let us return to the proof of the Wigner's law using the supersymmetric for- 
malism. The advantage of supersymmetry is to allow the integrated resolvent 
to be written as a functional integral in an ordinary field theory because the 
corresponding normalizing determinant is 1. More precisely let us return to 

< v{E) >= - hm f P(H)dH^-%Tr 1 . (2.32) 
e->o + J nN L — H + is 

We would have indeed with ordinary variables 

< (E-H + ie)^ >=-i J P(H) dH J )% E _ H+is)Sd <Zs (2 ' 33) 

where the imaginary sign has been chosen so that the bosonic part of the integral 
converges, using the small positive part e~ eS s . The integral of such a quotient 
is not easy to manipulate. But with supervariables taking the self-normalization 
into account we can write 

(E-H + te)^ 1 = -i J S*Sje^^ E - H+ie ^d^*d^ (2.34) 



G 



so that integration over H can be performed explicitly. Indeed 



< (E — H + le)^ 1 >= -i J P(H)dH J S*Sje l<s> ^ E - H+tE ^d^*d^ (2.35) 

and the Gaussian integral of an exponential linear in H can be performed ex- 
actly, giving rise to an exponential quartic in <I>, something which physicists call 
a " $ 4 interaction " : 

< ( ' N; "'" "' >= exp{-(l/2JV)^($?$ j )(*5*i)}- ( 2 - 36 ) 



The essential step is now to recast this quartic term in the form of a so-called 
"vector model", that is to factorize the quartic term as a square. Taking care- 
fully into account the anticommutation rules one finds: 

$,•)(*,• $0 = {Y.StSif + 2(£ S?Xi) - (ExiXi) 2 - (2-37) 



When rewritten this way, one can introduce a representation of the quartic 
term as an integral over a single supermatrix field R. Physicist call this idea 
a " Hubbard-Stratonovich" transformation. In practice the many variables H 
have been reduced to four "mean field" variables, the coefficients of R. We see 
that the key fact under this phenomenon is the independence of the variables 
Htj 1 . 

More precisely we have 
ex P{-2^E( $ *^)( $ 5*0}= [ dRexp{-jStrR 2 -iJ2$*R$i} (2.38) 

ij ' " i 

with 

Indeed due to the i factor in the fermion-fermion b term, the bosonic part 
of the quadratic form Stri? 2 = a 2 + b 2 + 2p* p is positively definite, ensuring 
convergence. 

We can now interchange the $ and R integrations and perform the $ inte- 



1 For a random matrix whose coefficients are not independent, this mean-field "Hubbard- 
Stratonovich" transformation is not valid and Wigner's law may no longer hold; however if 
a random matrix has for instance N 2 /p(N) independent coefficients, each spanning an orbit 
of p(N) sites in the matrix, we may search for a reduced representation with only 0(p(N)) 
"Hubbard-Stratonovich" variables and still expect Wigner's law for instance if p(N)/N — > 
as — ► oo, see below. 
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grations which are now quadratic: 
< Tr(E - H + ie))- 1 > = -i J dR J d$*d$^S* k S k 



exp{i(E + ie) - yStri? 2 - i ]T 

i i 

= N j dR(E — R + ie)^ 1 

exp{-yStri? 2 - NStv\og(E - R + ie)} , (2.40) 



so that 



< v{E) > = - lim ^3? J dR(E - R + ie)^ 1 

exp{-yStri? 2 -iVStrlog(S-i? + ie)} . (2.41) 

This R integral as N — > oo should be accessible to saddle point analysis, 
since the integrand is of the form e~ NS ( R ^ with N large in front of the action 
S(R) = iStri? 2 + Strlog(£ - R + ie). The saddle points satisfy in the limit 
e -» 

R= (E - R)- 1 : R 2 - ER+1 = (2.42) 
so the corresponding eigenvalues are 



£ = E/2- iy/i -E 2 /A , £* = £/2 + iy/l - E 2 /A, (2.43) 
with EE* = 1, E — £ = £* . The corresponding bosonic saddle points are 



(2.44) 



with s a = ±1, s;, = ±1 These saddle points are not along the original path of 
integration, which is a real and b real. So we have to shift the a integration 
contour to include an imaginary part —is a y/l — E 2 /4 and the b contour to 
include an imaginary part —iE/2. We remark that in order to avoid crossing 
singularities we must take into account the positive sign of e and choose s a = +1 
if e > 0. So we are left with two saddle points at Sb = ±1. But it is possible 
to explicitly perform the Fermionic integration and check that the saddle point 
at Sb = 1 dominates the saddle point at Sb = — 1, which is suppressed by a 1/N 
factor. 

Indeed we have 

"-^'-^"- ^-.Kg-a) '"'- (2 ' 45) 

SM{E - R) = §^ 6l l- l ,- fi (E _J (E _ ib) ] . (2.46) 
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Hence 



< v(E) > = I dadbdp*dpe' h < 1 



7T 



(E-ib)» M . 1 



so that, performing exactly the p and p* integration, and translating the con- 
tours, we get, remembering that by l|2.2I[l and anticommutation J dp*dpp*p = 



j_ 

2ir ' 



<u(E)> = -L$*/ da f dbe- N ^ a2+b ^ 



a-a=- x /l-£; 2 /4 JQb=-E/2 



(E ib) 1 + o(l/N)} . (2A8) 



(E-a) N + 1 L (E-a)(E-ib) 

At the vicinity of the first saddle point we have a — £ + 5a, ib = £ + iSb and by 
ordinary Hessian analysis we need to expand to second order in 5a and 5b. We 
find 



c -N/2(a 2 +b 2 ) Cg lb ) = e -N/2(l-£ 2 )(Sa 2 +6b 2 )+0(S*) _ ^ ^ 



N 



(E-a 



Furthermore at the first saddle point we have [— 1+ (E-a)(E~ib) \ = — ^ +£ 2 and 
( g _ a ) = 5 so that the Hessian approximation near the first saddle point gives 
a contribution 

< v x {E) > = ^ J 5a J 5be- N ^ 1 - £2 ^ Sa2+Sh ^£(-l + £ 2 + OQ./N)) 
= -^l-E 2 /4 + 0(l/N) . (2.50) 

7T 

At the second saddle point we have a = £ + 5a, ib = £* + i5b so that — f + 
(E-a)(E-ib) = a t this other saddle point. Therefore to leading order as N — ► oo 
there is no contribution to v(E) from this other saddle point. It is possible to 
bound the contributions away from the saddles and to conclude therefore 



<^>>= r™ 3 ,aiEi< - 2 ' (2.51, 

I otherwise . 

The computation of level correlations is similar but complicated by lack of 
positivity of the Hubbard-Stratonovich form, which leads to technical compli- 
cations. There is also an additional symmetry, and the main non-trivial mean 
field integral therefore takes values in a coset superspace parameterized by eight 
real variables, four bosonic and four Fermionic 6 . 

Let us also stress that the corrections as N — ► oo to Wigner's law can be 
in principle systematically computed but the computation becomes more and 
more difficult in practice for higher order terms. However it is relatively easy 
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to derive e.g. crude bounds on the probability that a random N x N matrix 
develops a large norm using some kind of Tchebychcff inequalities. For instance 
it is an instructive exercise to prove |5] 

Lemma 1 For sufficiently large a and N one has 

Prob(||iJ|| > aV6) < 8Na 6 /27e- Nl/3a2 /3 . (2.52) 
The Tchebycheff inequality here is simply to use 

Prob(||iJ|| > aV6) < {6a 2 r m < TrH 2m > (2.53) 

and to optimize over m. Gaussian integration for < TrH 2m > can be explicitly 
performed through Wick's theorem. It leads to the evaluation of so called Feyn- 
man graphs. A careful counting of these graphs leads to the proof of the lemma 
by taking m proportional to a 2 . This bound is not optimal but shows clearly 
that the probability for a random matrix normalized in this way to develop 
eigenvalues much larger than O(l) is very small. Such crude bounds are suffi- 
cient for complete control of the Anderson model in the regime \p 2 — E\ » A 2 , 
see below. 

3 The Anderson Model 

The Anderson model of an electron in a random potential corresponds to the 
Hamiltonian 

H = -A + XV(x), 

acting on the Hilbert space L 2 (R 2 ), where A is the usual Laplacian and V is 
a real Gaussian process on R 2 with short range correlations. A is a coupling 
constant that allows to adjust the strength of the random disorder with respect 
to the deterministic part. 

This model was initially proposed to describe the electron motion in doped 
semi-conductors at low temperature or in normal disordered metals. It is now 
the central model for the theory of electronic transport and wave propagation 
in disordered systems It was conjectured by Anderson as soon as 1958 [TU] 
that such a model exhibits a localized phase in which the electrons are trapped 
by the defects. In 1979 it was argued that this model has a phase transition in 
dimensions three or more between the localized phase and an extended one jll| . 

The localized phase is now well under control. In one dimension, localization 
was rigorously established for any disorder at the end of the seventies |12l 113) . 
Later localization was established in any dimension at strong disorder or for 
energies out of the conduction bands ^3] . A simplified and more efficient method 
to get this result was given in |15j . 

In contrast the weak disorder regime is still poorly understood. In two di- 
mensions it has been argued |1 II I16| and numerically established 17 that local- 
ization persists at arbitrarily small disorder, with a localization length diverging 
like 0(e c / A ). In dimension three, numerical simulations confirm the existence 
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of the Anderson transition |17|. leading to an extended phase. In addition, ana- 
lytical results ^Hl an d other numerical calculations show that the level spacing 
distribution follows the Wigner-Dyson distribution for random matrix theory 
(RMT) [PQ. This gave the motivation for a description of mesoscopic system 
in terms of RMT [201 - This method has been very successful when compared 
to experiments and was the source of developments of supersymmetric meth- 
ods in solid state physics. But this heuristic connection between random 
matrices and the Anderson model remained mysterious. And on the rigorous 
level it is still a mathematical challenge up to now to prove even regularity 
of the DOS at weak disorder in the conduction band (that is, analyticity in 
energy in a narrow band around a real interval in the conduction band of the 
deterministic Hamiltonian, see Conjecture 1 below). This is still unknown in 
either dimensions two or three (see |22| for recent continuity results outside the 
conduction band). 

Using a phase-space analysis inspired by the renormalization group method 
around Fermi surfaces in condensed matter G. Poirot and coau- 

thors j2H] have understood the connection between the Anderson model and 
random matrices. They established that the effective Hamiltonian near the 
Fermi level is given indeed by a random matrix model. In the simplest case, 
that of two dimensions this random matrix model is almost identical to the GUE 
(Gaussian Unitary Ensemble), but just contains one extra discrete symmetry. 
This symmetry is called the Flip Symmetry in pp, and the corresponding Flip 
Matrix model has been analyzed through the supersymmetric method. In three 
dimensions the flip symmetry becomes a continuous U(l) symmetry, and pro- 
duces more complicated correlations between matrix elements |2(i| ; however even 
in that case, these correlations are not expected to change the main statistical 
properties of the spectrum, such as Wigner's law for the density of states. 

In the coming subsections we will summarize the content of P^. Let us 
recall that in the flip symmetry was slightly simplified. This flip symmetry is 
indeed essentially a Z 2 symmetry, but which effectively degenerates into a larger 
symmetry near the diagonals of the matrix, which correspond to degenerate 
rhombuses in the momentum representation. Therefore in the last section of this 
review we will refine the phase space analysis of the Anderson model and propose 
improved matrix models to study this slight degeneracy near the diagonals of 
the flip symmetry, using some hierarchical approximations that mimic the true 
situation but are more tractable from the point of view of the supersymmetric 
mean field analysis. Similar ideas could also apply to the study of the U(l) 
three dimensional symmetry. 

3.1 Phase Space Analysis 

As in other condensed matter models, the ultraviolet region in momentum space 
is irrelevant and cut off in the Anderson model. The most common way to do this 
is to replace ordinary space by a lattice such as Z d . This corresponds physically 
to the so-called tight binding approximation in which the electron spends all its 
time on lattice sites in a crystal, jumping from one site to an other. This lattice 
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model automatically cuts off large momenta since the Fourier transform of the 
lattice Z d is simply the torus [— 7r,7r] d . 

Since the lattice cutoff breaks rotation invariance, it is also common to con- 
sider a rotation invariant momentum cut off. For that we consider the free 
Hamiltonian without the random potential, which is simply the Laplacian in a 
suitable unit system: 

H = -A (3.54) 
and the free retarded Green's function: 

Go = Ho-E-ie - (3 ' 55) 

also called the propagator in physics. In Fourier space this is a diagonal multi- 
plication operator. A rotation invariant momentum cutoff can be implemented 
by multiplying this operator by a smooth C°° function k with compact support 
which suppresses large momenta. For instance we can ask that this function is 
identically 1 for \p\ < A and for \p\ > A + 1, where A is a fixed number, e.g. 
10. This leads to define the cutoff free retarded Green's function: 



«(p) 
p 2 — E — is 



GoAp) = ~ 2 V - ■ (3.56) 



This is usually called in physics the "jcllium" cutoff, and we call the corre- 
sponding model the jellium Anderson model. Regular perturbation theory for 
this model consists in computing the interacting Green's function as a power 
series in A through a resolvent expansion: 

G = Go , : * = Go - GoWGo + G XVG XVG + ... (3.57) 
1 + A V Go 

The rotation invariant cutoff consists in replacing everywhere Go by Go, K . Re- 
mark that the convergence of the power series for fixed V is no problem at A 
small enough, since this is a geometric series. However we are not interested 
in any particular V but on the average over V. The corresponding integrated 
Green's function is obtained by integrating over V, using the Gaussian rules of 
integrations, called Wick's theorem in physics. Only terms with even numbers 
of V survive and the result can be indexed by drawings called Feynman graphs. 
In these Feynman graphs, the Go propagators can be pictured as full lines, 
whether the covariances for contracted V, which are S functions in x-space, are 
usually represented by dotted lines (see Figure 1). In this case the first non 
trivial graph also called tadpole graph (see Figure 2) gives a momentum space 
contribution 

G = A 2 [d d p 2 K{ 1 ] , . (3.58) 
J p — E — l£ 

This is a pure number, with a finite imaginary part 5 as e — > 0. For instance 
in two dimensions and for E = 1, this imaginary part S is A 2 7r 2 . Therefore 
physicists, resumming only the graphs G n with n consecutive tadpoles (see 
Figure 3), which form a geometric series, and neglecting all other graphs, can 
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Figure 1: A typical Feynman graph 



Figure 2: The first non trivial graph, the "tadpole' 



perform the limit e — > 0: 



DC 



k(p) 



G tadpole — Go / y (5Gp) 
n=0 



n 



V 



i 2 — E — iSn(p) ' 



(3.59) 



With their usual boldness, they conclude that the averaged Green's function 
should be regular in momentum space, and consequently also decay in x-space 
at a spatial scale inverse of 8, at least for small A. The corresponding conjectures, 
alas unproved yet, could be described as follows: 

Conjecture 1: Analyticity of the Averaged Density of States 

Let [a, b] is a fixed interval well inside the spectrum of the unperturbed Lapla- 
cian spectrum (this means [a,b] c]0, 2d[ for the lattice model, and [a,b] C]0,A[ 
for the jellium model). For any X small, the averaged density of states < v(E) > 
is analytic in E in a rectangle centered on the real interval [a, b] and of width 
0(X 2 ), hence of type [a,b] x [-iO(\ 2 ), +iO(\ 2 )]. 

Conjecture 2: Scaled Polynomial Decay of the Averaged Green's 
Function 

There exists a constant c such that for any integer q, there exists a constant 
K q such that the averaged Green's function < G > in direct space obeys long 
range spatial decay of power q and scale 0{\ 2 ): 



These conjectures arc not completely optimal, but even in this form they are 
still a challenge. It is in fact expected that we have analyticity in a rectangle 



\<G(x,y)>\ < K q (l + c\ 2 \x - y\)-« . 



(3.60) 



Figure 3: The graph G n with n tadpoles 
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with [a, b] closer and closer to the full unperturbed spectrum as A — * 0. Also 
in the case of the lattice model one expects in fact exponential scaled decay 
of order A 2 . Since a compact support cutoff k cannot be analytic, its Fourier 
transform does not decay exponentially. So for the jellium model one expects 
at most fractional exponential decay (if the function n is Gevrey) . 

Perhaps more interesting than the conjectures themselves is what they mean: 
namely that no matter how small the interaction, it regularizes the average 
Green's function. This is somewhat similar the phenomenon of dynamic mass 
generation in field theories like QCD, although here the regulator is polynomial 
in A. 

These conjectures would be important also to justify perturbation theory and 
the tadpole approximation in the small coupling regime of the Anderson model. 
As is well known the ordinary perturbative approach is plagued by mathematical 
problems. Indeed due to the large number of Feynman graphs (there are nil = 
(n — l)(n — 3). ..5. 3.1 graphs at order n) the integrated perturbation series has 
zero radius of convergence. Over the years, techniques have been developed to 
overcome this problem by replacing the complete divergent Feynman series by 
carefully truncated expansions with Taylor remainders that can be rigorously 
bounded. This set of techniques goes under the name of "constructive field 
theory" |27) . 

Constructive field theory tells us that in order to achieve non-perturbative 
theorems in this kind of situations where the interplay between propagation 
and interaction is non-trivial, the best frame is to analyze the theory in phase- 
space. Some multi-scale renormalization group analysis is usually necessary. 
Renormalization group is the tool which computes long range behavior from 
local physical interactions. Cluster and Mayer expansions are the correct non- 
perturbative steps which make the renormalization group mathematically well- 
defined [23 - The guiding line is to identify the relevant degrees of freedom 
of the problem in phase space, and then to perform some carefully controlled 
perturbation theory which roughly derives at most one perturbation step for 
each such degree of freedom. 

A Hamiltonian, whether random or not, acts indeed on a Hilbert space of 
quantum states. Using a system of units in which h = 1 a quantum degree 
of freedom corresponds roughly to a unit volume in phase space. Although 
the quantum Hilbert space for a particular problem may be formally infinite 
dimensional (such as L 2 (R)...), in practice any physical problem should in fact 
be well approximated by a finite dimensional version of this Hilbert space. The 
main problem is to identify the relevant degrees of freedom of the system, and to 
control their proliferation in the ultraviolet (field theory) or infrared (statistical 
mechanics) limit. 

A partition of the classical phase space into cells of unit volume may be 
viewed as giving a particular basis of this finite dimensional Hilbert space. Not 
all partitions are admissible; they should roughly respect the Heisenberg uncer- 
tainty principle ultimately related to the symplectic geometry of phase space. 
Any product of the typical dimensions of any cell in the directions of a conjugate 
pair of positions and momenta should not be much smaller than 1. 



14 



Some decompositions may be more convenient than others for particular 
problems. In the lattice Anderson model on a lattice, the phase space is 1 d x 
[0,27r] d since positions lie on a square lattice and momenta lie in the dual tori. 
In the jellium Anderson model, it is M. d x B(0,A + 1), where the closed ball 
B(0,A + 1) in K d with center and radius A + 1 is the compact support of 
the cutoff function k. This is not very different from the previous case since k 
restricts the large momenta, and the x-space should be again well approximated 
by a lattice of cubes of side size roughly the inverse of A, hence of order O(l). 
Indeed all propagators and functions built out of them will be roughly constant 
within each of these cubes. 

The initial natural partitioning of such a phase space is to associate to each 
lattice site or cube in direct a;-space a degree of freedom. Momentum space 
is compact, and the only reason for which the problem has infinite number of 
degrees of freedom is the infrared problem, or thermodynamic limit: we want 
to understand the theory in infinite volume in a;-space. 

This obvious partitioning is well adapted to the study of the localized phase 
of the Anderson model, in which the random potential dominates over the prop- 
agation. Indeed a random identically distributed potential V is simply a mul- 
tiplication operator in x-space. It is a diagonal operator in x space with i.i.d. 
diagonal coefficients V(x) for x € Z rf . Spectrum is pure point and the eigen- 
values of this operator are obviously the wave functions concentrated on single 
sites, or 5 functions in x space. The density of states and statistics of eigenvalues 
are Poissonian: the presence of an eigenvalue at some energy E simply means 
that for some x we have V{x) = E, and this information has no influence on 
the values of V(y) for y ^ x. Therefore in this case i?2(s) ~ s ~o S(s) + 1. 

Most of the rigorous work on the Anderson model is concerned with the 
localized phase. In this phase the addition of a hopping term such as a lat- 
tice Laplacian is a small perturbation which does not fundamentally modify the 
spectrum, which remains pure point, and the localized character of the eigen- 
f unctions. 

But this partitioning of phase space is not the right one to study the other 
regime in which it is the lattice Laplacian with its continuous spectrum which 
is the main effect, and the random potential which should be the perturbation. 

From now on let us restrict to a single averaged Green's function or the 
density of states, leaving the level correlations for future studies. To prove 
analyticity in the rectangle of Conjecture 1, we need only to consider a fixed 
value Eq in [a, b] and, translating E to Eq + E, to prove analyticity for \E\ < 
0(X 2 ). If we can do this for any E e [a, b], Conjecture 1 follows. Since this 
will not change anything except a global rescaling we always choose Eq = 1 in 
what follows. 

If we take seriously the propagator as a guide to the right partitioning of 
phase space, we recognize that the interesting region in momentum space is 
the singular region p 2 = 1. This is the d-dimensional sphere, for instance 
in two dimensions it is the circle of radius 1. We know from the theory of 
rcnormalization group around extended singularities |241 |2*B] that in this case 
we should perform a change of basis, dividing the region p 2 ~ 1 into smaller and 
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smaller cells to analyze it with great care; correspondingly x space should be 
cut into dual boxes that will be larger and larger as we approach the singularity, 
hence correspond to longer and longer distance effects. 

This decomposition correspond to cut first momentum space into shells ac- 
cording to a geometric progression which pinches more an more the singularity. 
Fixing a rate M > 1 for this geometric progression (for instance M = 2), this 
means that we write 

«(p)=5>i(p) (3.61) 

3=0 

where Kj (p) is a smooth function with compact support which roughly ensures 
\p 2 - 1| ~ M~i . 

The index j plays the role of a renormalization group index, and the cor- 
responding region in momentum space is called the j-th slice or momentum 
shell. The last index j ma x is chosen so that M _ -' max ~ 5 = 7r 2 A 2 . Indeed the 
expansion should not be infinite. The last function Kj max (p) for the last slice 
jmax should be different, meaning only \p 2 — 1| < M^ ima " . The reason is that 
we expect the basic physical intuition behind 1)3. 59[) to be correct. Averaged 
Green's functions should decay after that scale; equivalently the singularity of 
the free propagator should be screened by the appearance of an imaginary part 
due to the interaction with the random potential. 

This screening phenomenon, like mass generation under influence of interac- 
tion in field theory, can only occur at a scale where interaction effects become 
of the same size than the free propagator. The conclusion is that we must 
distinguish two regimes: 

- The regime \p 2 — 1| >> A 2 , or j < j max . In this region the free propagator 
should dominate over the interaction and regular cluster and Mayer expansion 
should prove that this part of the Green's function remains in every respect 
close to the free one. 

- The regime \p 2 — 1| >> A 2 or j ~ j max , which corresponds to a few 
renormalization group slices (or in fact just one if we take M big enough). In 
this last slice clearly ordinary perturbation theory (even recast under the form of 
cluster or Mayer expansion) cannot work since propagators and typical random 
potentials are of the same size (otherwise the non-trivial effect of screening 
would be impossible). 

In this regime the functional integral oscillates wildly. A first attempt to 
exploit the corresponding cancellations through Ward identities was made in 
30 . However this turned out to be a difficult and not completely conclusive 
approach. Oscillating integrals with analyticity properties can often be better 
analyzed by finding saddle points and shifting contours to pass through them. 
In favorable cases the leading part of the integral then comes simply from a 
Hessian approximation near a single dominant saddle. This saddle becomes a 
new expansion point for the theory, similar to a change of vacuum in field theory. 
The corresponding shift is definitely a non-perturbative tool, very difficult to 
understand in terms of combinations of the initial Feynman graphs of Figure 1. 
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Such a shift precisely occurs in the supersymmetric formalism of the previous 
section: the contour for the mean fields a and b has been shifted to pass through 
a non trivial saddle point that lied out of the original integration contour. In 
this way a highly oscillating integral, has been replaced by an equal one which 
is much simpler to analyze. We must explain now why the Anderson model is 
analogous to a matrix model, so that the same strategy should also work there. 

3.2 Sectors and the Random Matrix Approximation 

Momentum slices or shells without further decomposition are not suited for 
a phase space analysis. Indeed a free propagator restricted to such a shell 
has no simple dual decay properties in direct space. The correct phase space 
analysis must split further each shell into angular sectors, along the directions 
tangential to the extended singularity. What limits the length of the sectors 
in these directions is the curvature of the singularity (see also |2H] for a 
non-spherical example). In the case of a spherical singularity such as the one of 
the jellium Anderson model, there are two natural and acceptable solutions to 
the splitting into sectors: 

- split the spherical shell of width M - - 7 into roughly M^ -1 ^ sectors (noted 
by Greek letters such as a or a) with same size in all directions. This is the most 
natural idea and the corresponding sectors are called isotropic. It is however 
not optimal for momentum conservation analysis. The direct space is then cut 
into a lattice T>j of cubic boxes A of side size M- 7 for each value of j. 

- split the spherical shell of width M - - 7 into roughly M^ -1 ^/ 2 sectors, still 
noted a, of length M~i/ 2 in all directions tangent to the singularity. These 
sectors are called anisotropic. This is the optimal tool to analyze e.g. momentum 
conservation, but the price to pay is that the dual decomposition of x-space gives 
rise to anisotropic parallelepipeds which change as the sector changes. For each 
sector a in a shell j we have a different lattice T> a of spatial boxes covering 
R d , which have length AP in the direction of a, and length only AP' 2 in the 
other orthogonal directions (since for the sphere these orthogonal directions 
correspond to the tangent directions at a). 

In the simple isotropic case, a phase space cell of roughly unit volume is 
therefore made of a value of j in [0, j maa: (A)], a spatial box A S T>j and a 
momentum cell a in the j-th shell. The true Hilbert space of the theory is well 
approximated by a finite dimensional Hilbert space whose basis is indexed by 
these phase space cells. 

Now at last comes the analogy with random matrices. In x space, V is a 
real diagonal operator; therefore in p space V(p) is complex with 

V(-p) = V{p), (3.62) 

and V is a random convolution operator Vf(p) — J dqV(q)f(p + q) with covari- 
ance < V(p)V(q) >= S(p — q). A delta or short range covariance in x space 
translates into a white noise in p space, that is all momenta correspond to in- 
dependent variables with a priori equal weight. If we use isotropic sectors to 
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simplify, it is natural to discretize the momentum space for V with cells of the 
same dimensions that the angular sectors used for the propagator. When two 
propagators with momenta a and j3 join at a potential V, the momentum of 
V, or "transfer momentum" by overall translation invariance of our theory, is 
p = a — p. There is therefore an effective ultraviolet cutoff on V (because a 
and (3 being bounded, so is a — (3) . 

Here is a first subtlety. There is no reason for a and to lie in the same 
shell. The operator V can be decomposed into the shell-conserving or diagonal 
part V s — V s ,j and the shell-changing or off-diagonal part V . For the time 
being, let us neglect V and consider an operator V s j acting in a single shell. 

The covariance for V has no singularity and we should therefore discretize 
the momenta for V into roughly M dj cells r of side size . Hence if we divide 
space according to the lattice Vj , to a given cube A should correspond roughly 
M"i discrete variables V t .a, which are roughly independent and identically dis- 
tributed. The interaction being local in x space, such a variable V t .a must join 
sectors (a, A) and (/3, A) spatially localized in the same cube A of T>j. From 
now on fixing such a cube A, we can therefore consider the variables V T ,A, as 
the discrete elements of a random matrix that we call V(A). 

Let a and (3 be two indices parametrising two sectors of the same shell j. 
We need to know which pairs of sectors a given independent variable V T) a can 
join, so to analyze the relation r = a — j3. The answer happens to be strongly 
dimension-dependent . 

A first remark is that if r = a — (3, then — r = (3 — a. This means, using the 
relation (|3.62|l that the matrix V must be hermitian, in any dimension: 

V a ,p = Vp, a . (3.63) 

In two dimensions if we consider a given r of length between and 2E, there are 
typically roughly two pairs of momenta of length E, forming a rhombus such 
that t — a — (3. These two pairs are (a, (3) and (—/3, —a). Therefore the matrix 

V in addition to being hermitian has an additional symmetry 

V a ,f3 = V- ,- a (3.64) 

called the flip symmetry in *Q- Under this symmetry, two coefficients in e.g. 
the upper triangle of the matrix correspond to the same random variable. This 
means that V is not exactly a GUE. We can label the sectors of the circle in 
such a way that these two identified coefficients are symmetric with respect 
to the antidiagonal of the matrix. For that it suffices to divide the circle into 
an even number 2N of sectors labeled as {1, • • • , N} U {— N, • • • , —1}, so that 
{1, • • • , N} labels projective sectors, as shown in Figure 4. Of course this is only 
a generic situation, and when the rhombus degenerates the answers are different. 
For instance for a momentum r of length 2, there is only one pair (a, (3) = 
(t/2, — t/2) such that r = a — (3, and for a momentum r of length 0, every pair 
(a, a) is a solution. This means that all coefficients on the diagonal of the matrix 

V are equal. Intermediate situations when the rhombus almost degenerates 
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Figure 4: Numbering sectors 



should also complicate the flip symmetry, and these issues are discussed in some 
detail below. 

In three dimensions for a given r there is typically an orbit of y/~N pairs 
of sectors (a, (3) such that r = a — /?, forming a cone. The discrete Z 2 flip 
symmetry is therefore replaced by a U(l) symmetry; the set of upper diagonal 
coefficients of the matrix V are therefore split into 0{N Z / 2 ) orbits under this 
symmetry, each orbit containing roughly >/~N coefficients. Clearly this matrix 
model is farther from the GUE model, but still the density of states should fall 
in the same category. 

Having recognized the analogy, in phase space, between the Anderson model 
and random matrices, we shall now review first the regime \p 2 — E\ » \ 2 ir 2 , 
then the last shell \p 2 - 1| ~ \ 2 ir 2 . 

3.3 The regime \p 2 - 1\ » \ 2 n 2 

This is the region in phase space where by definition the deterministic part — A — 
1 is forced to be not too small, and where the random part XV is statistically 
very small compared to this deterministic part. However the distribution for V 
is not compactly supported. Therefore the probability for V to have a norm 
comparable in size to —A — 1 is never zero, even in that regime. However it is 
so small that if such an event happens in a phase space region R, there is such 
a small associated probability factor that it can pay for a very crude bound, 
obtained by translating directly the integration contour for the potential V in 
that region R. In this way we can prove not only that such an event has small 
probability, but that it does not influence other regions of phase space. 
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The corresponding rigorous analysis has been performed in great detail in 
the work of G. Poirot [2H , and here we give only a rough sketch of the argument. 

To get an idea of the respective sizes of the deterministic and random pieces 
in (H — 1) we can first work at fixed j for a shell-conserving piece V Sj j. Consider 
first the deterministic part p 2 — 1 . Restricted to the j-th shell it is by definition of 
that shell an operator with a norm of order M~ J . We want to prove that when 
restricted to that shell, the random part has typically a much smaller norm. 
If this is the case ordinary perturbation theory of the random part should be 
trusted. 

The full interacting averaged Green's function is 

s , <p l . (i + 2 . r . (a.es) 

p 2 - 1 - E + is p 2 - 1 - E + is ' v ; 

We have divided into shells through k = J^j K j ■ F° r a slightly more symmetric 
form we can put square roots of the free propagator on each side of V, and write 
Kj = j] 2 . Then the random operator restricted to the j'-th shell is 



Vj(p) ^ Vj(l) 

(p 2 - 1 - E + ie) 



,-"2.^ V<,-i) ™ (3.66) 



with 

V(p-q) = I dxe lx -( p -^V(x) . (3.67) 



The main problem to define and control the averaged Green's function is now 
to compare the norm of Hj to 1 in order to invert the operator 1 + Hj . 

The sectors in a cube are roughly L 2 normalized characteristic functions 
e a — M d ^ 2 Xa{p), where the factor M d ^ 2 is required for normalization. 

The operator Hj restricted to a spatial cube of size Mi is now analogous to 
a matrix between sectors with coefficients 

H a>/3 ~ <e a ,He p > (3.68) 

Xaip) , T// v xp(q) 



2 



(p 2 - 1 - E + ie) 1 / 2 ^ H ' {q 2 - 1 -E + ie) 1 / 

since rjjXa — Xa, VjXis — Xp- I n a shell of index j, we have roughly (p 2 — 1 + 
ie) -1 / 2 = Mi I 2 , so that the covariance between coefficients can be computed as 



<H a ^,H^ s>v ~ M( 2d+2 WA 2 / dpdqdp'dq' Xa {p) X[ Mx^p')xs{q') 

dxdye ix ^ p - q h- iy -^'-^ < V(x)V(y) > v . (3.69) 
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Using < V(x)V(y) >v= 5{x — y) we find 

<H a ^,H J ^ s > v ~ M {2d+2)j \ 2 J dpdqdp'dq 1 

Xa(p)Xl3{g)Xj(p')Xs{q')S({p - l) - (p' - <?')) 
„ M (2d+2) 3X 2 I d P dqd P ' Xa (p)Xfs(q)x 7 (p')X6(p'-p + q) 
~ A a (<y a , 7 */J,« + *a,-«*/J,- 7 ) ford = 2. (3.70) 

In the last line we specialize to the generic momentum conservation rules for 
d = 2 (non-degenerate rhombuses) 2 . 

The random operator Hj in two dimensions restricted to a single cube, is 
therefore similar to a random matrix of the GUE type with the additional flip or 
U(l) symmetry, and with covariance proportional to A 2 . The size of the matrix 
is N ~ AP , the number of isotropic sectors. We know the average norm of a N 
by N GUE matrix to be not much more than 0(1) for matrices with covariance 
1/N on the independent coefficients (see Lemma 1 in section 2). Rescaling to 
covariance A 2 is like rescaling H by Xy/W ~ AM 3 ' 2 . Hence we conclude that 
statistically the norm of the random part Hj should not be much larger than 
0(1)AM J '/ 2 . This means it is statistically much smaller than 1 for AM 3 / 2 « 1, 
which is nothing but j < j ma x(A) (if M is large). We conclude that statistically 
we can perturb the random part with respect to the deterministic part in the 
whole regime \p 2 — 1| >> A 2 , as announced. 

Of course this is a very rough argument and we must provide details on what 
to do in the infinite volume case, and how to treat the statistically rare events 
where the random part is unfortunately much larger than expected. 

Since the space volume can be large, one performs a battery of tests which for 
each cube A tell us whether the random part in that cube A has an exceptionally 
large norm or not. This battery of tests is called a large/small field analysis in 
constructive theory. 

Roughly speaking one uses then an analyticity argument in the cubes where 
V has anomalously large norm, which form the so-called large field region. The 
contour of integration for the V's of this large field regions is translated to create 
an imaginary part of same order than the deterministic part, and a bound 
is applied to a resolvent expansion which tests whether the Green's function 
effectively or not has visited these large field regions |29| . 

Taking absolute values for this bound destroys oscillations so we need to 
pay for a large factor corresponding to the apparent change in the Gaussian 
normalization: 

\ e -(V +t af/cy e -V*/c = & W , >c (3 n) 

but this large factor is compensated by the small probabilistic factor coming 
from the large field condition. To understand roughly why this is possible in 

2 For general d we would find a factor A 2 M( 2-ti ) J times the constraint that there are 
momenta p, q and p' in the sectors a, (3, 7 such that p' — p + q is in <5, and this is a more 
complicated constraint. 
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two dimensions, where N ~ M J , recall that to create with XV an imaginary 
part of size M~ J requires a translation of the zero momentum Vq coefficient on 
the diagonal of the matrix V by iA _1 M _J . Then the factor to pay following 
<|3.71fl is e a l c with a = A _1 M - ^ and the covariance c for the zero momentum 
coefficient is M -2 - 7 . So the factor to pay is e A 2 . But if the norm of V is larger 
than A^M - - 7 ', by lemma 1, the small probabilistic factor is of order e~ M3 x ^ 
so it compensates for much more than the factor to pay. This is what essentially 
fuels the analysis of 29 . 

This large field analysis is also complicated by some auxiliary expansions 
to take into account the non-diagonal shell changing parts of the potential, 
which does not change the overall picture. Altogether this analysis, although 
technically tedious, is obviously much cruder than a shift of the expansion point, 
which becomes necessary to treat the second regime. 

The conclusion is that in this regime one can invert the operator H — E, and 
the averaged Green's function restricted to that regime has the expected size 
and spatial scaled decay rate in \xM~^ mam \ = \ 2 \x\. 



4 The regime \p 2 — 1| < cA 2 

We know heuristically from first order perturbation theory that the retarded or 
advanced Green's functions should decay with a rate of order A 2 7r 2 . Therefore 
it makes sense, at the level of a single Green's function or for the density of 
states to treat this regime as a single momentum slice. Expecting this decay, 
it is enough to understand the model in a single space cube of side size A -2 , 
and then to treat the full model through some cluster expansion. In this section 
we want to summarize the more detailed computations in £Q which prove that 
the density of state of a single matrix with the I2 flip symmetry has the same 
large N limit than a GUE matrix. The core of the argument reproduces the 
supersymmetric computation of section 2.3, but with some additional auxiliary 
fields to take into account the flip symmetry, ft is then shown that the presence 
of these fields do not modify the main saddle point and the leading contribution 
to the density of states. 



4.1 The Flip Matrix model 

The main idea in pQ is to introduce a discrete version of this single cube prob- 
lem, in which phase space and the 2d momentum conservation rules have been 
discretized and somewhat simplified, according to the isotropic discretization 
and to equations (|3.63I3.64|I of subsection 3.2. 

Using the labeling of Fig f , the structure of V can be summarized as follows 
(remember rows and columns are labeled as {1, • • • ,N, —N, • • ■ , — 1 } hence not 
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as {!,••• ,N,-l,--- ,-N}): 



V = 



v 



v 



«,-£> 



v 



and the main result of 1 is: 



Vo 



lb 



»,-/3 



I'd 



v 



v J / 



Theorem 1 /n the large N -limit, the DOS of this flip matrix model converges 
to Wigner's semi-circle distribution. The corrections to the limit are uniformly 
bounded as 0(1/ N) as N — > oo. 

Let us summarize now the proof of this theorem. 



4.2 The isotropic 2d flip model: Fermionic terms 

Introducing supersymmetric fields 

for a = ±1 and Greek variables such as a, (3 now running only within [1, .... N], 
we find the following formula for the density of states 

u{E) = lim —3 / -i- £ S+ a S aa e lL <> J] d^ a #™ d^V), (4.72) 

OB— >0+ 7T / ZiV ' — ' 

Q.cr Q.cr 

where the supersymmetric action decomposes as 

Lq = (Ao +B Q + C* ), 
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with 

A °= S { A ^,-/3(^^/3 + *i ff/ 3*-«) + ^- c -}> 
a</3,CT 

5 = ^Ay a ,_ a (*+*_ a + /i.c), 

a 

Co = (£ + AF )^^*aa ■ 

The Gaussian integration over the V variables will be performed except for 
Vq, leading to the following quartic action 

v{E)= Um -Im ( S+S a e- L X\<M+ a dH aa dix Q {Vv), 

ImE—>0+ 7T J AJ - 

Q.fJ 

L = \ 2 {A + B + C), 

A= J2 (*^/J + *i^*-a)(^*a + *ia*-a/9), (4-73) 
a</3,a 

a 

C = (-iA^E - *A- V ) ^ 

CK,cr 

Some tedious rewriting using the commutation and anticommutation rules even- 
tually leads to L = \ 2 (A + B + C) with 

A = (1/2) + ^t a ^-a)(^c + *±Jf-*fl), 

a 

B a = -2S+S- a Sl a S a - \ ^(*+ a * CTQ ) 2 




a a 



C = C= (-i\- 2 E - iX^Vo) (4.74) 

The term B being a sum of diagonal quartic terms, will be neglected for a while. 
For indeed it cannot be written as the square of a sum over a. However, the 
bosonic part of this term has the wrong sign, which may create some difficulties 
when performing the integration in (|4.72|) . This problem will be addressed later. 

The other terms can be reorganized so as to get the square of a sum over 
the a's by pushing the terms with index a on the left and the ones with index (3 
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on the right. This must be done with care according to the commutation rules 
for Bosons and for fcrmions. The calculation is tedious but straightforward and 
gives: 

A = (1/2) ( Y, SLS aa ) 2 - (1/2) ( xLx.a) (h.c.) 

a, a a, a 

+ 2 ( £ S+Sl a ) (h.C.) + ( J2 SaaXta) (h.C.) 

+ (Y s ^xt aa )(h.c). (4.75) 

The squares can be unfolded by mean of an integration over auxiliary gaussian 
fields. This amount to introduce two real fields ao and bo, one complex field a, a 
and two pairs of Fermionic fields £, £ and r], fj with Gaussian measure 



d(i,(ao,bo,a,a,€,£,r),r]) = e 



-VZ /2~al/2-bl/2-\a\ 2 /2-Cli-rfv 



dV da db d 2 a t 
x • • • x — tt-^ a£ d£ drj drj. 



Therefore 



2(2tt) 1 /2 

2 ^+ c > = J dp cxp (zA]T <f+i?$ Q ) 



(4.76) 



with 



Y = (ao + V + \~ 1 E) ]T S+ a S™ 

a a. a 

+ (ib + V + \~ 1 E) xLXaa + (aJ2 S <* S -* + h - C - 

+ (V J2 S ™X* a + h.c^j + (rf S aaX -a a + h.c)) , 



here $ is the supcrficld $+ = (S*+, S- a , xt > 2 X-a) an d i? is a 4 x 4 supermatrix. 
Setting A) = (ao + Vfc + A" 1 ^) and zB = (ib +Vo + X~ 1 E) we obtain the 
following representation 



v(E) = -^S J dulXR}^ [Sdet(i\R) 
The matrix R is a so-called 4x4 super-matrix 



-JV 



R = 



f A a -Z77 \ 

a A rf 

£ 77 zBo 

\ -it]* «B y 



A p* 
p B 



(4.77) 



(4.78) 



and the coefficient [Ai?] 6b 1 11 is a particular (bosonic) coefficient of the inverse 
matrix, the one corresponding to first row and column. In this computation 
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we use the fact that this coefficient is also equal to [Ai?] b6 1 22 and that we get 
therefore 27V equal terms from the sum over a and a in 14. 7211 . and cancel this 
sum with the 1/2N normalization factor. 

We perform a rescaling of all the variables by A -1 = y/2N and write 



v[E) 



-1^ 
7T ^20r 



daodbodadadVodt;* dt; drf dr\ 



-N(a 2 +b 2 + \a\ 2 + V^+2(ei+V , 'v)) 



[Ac 



1 



(re+^*»7)](iBo 



,2A< 



(4.79) 



defining ^0 = 00 + ^0 + E, iB = ib +V +E. 
We want to compute the fermionic integration 



dCd^drj*drj e - 2Ar («*«+"» 



(4.80) 



Note that 



-2iV(4*«+jj*77) 



(4.81) 



[4> 



1 

\a\ 2 



2 2 
iB iBq 



-N-l 



jV + 1 



{iB y 



-N-l 



(4.82) 



|2\JV+1 



An - lal 2 



(*Bo) 



(N+l){N + 2) 



Al- 



(iB ) 



All terms vanish through fermionic integration except their terms in £* £77*77. 
So H4.8()(l can be computed as 



If 



1 



^(Ag-H 2 ; 



7V+1 



JV+1 4JV + 2 8N{N+ 1) 4„ +4 ^ 2 



A 2 - |a| 2 (iJ3 ) 2 ^ - |a| 2 iS 



JV + 1 4A 



i-Ao-Bo \A%- l°l 2 iB o 



AN 



\2\N+l 



JV+1 (AN -2 8NA 



Al -14 



(iB ) 



tB 



AN 2 



(4.83) 
4iV 



iA B 
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V (E) = —3^5 / da db dadddVoe- N(a « +b « + ^' 2+v o\iB ) 2N A Q I f , (4.84) 



HE) = / d^(A 2 -\a\ 2 )- N - 2 (iB ) 2N - 2 (4.85) 



We have therefore 

7T ^2^7? 

hence 

4tt 7 /2 

[(JV + 1) ((27V - 1)Aq - UNAqBq) + 2iNB (A 2 - \a\ 2 )(l + iNA a B a )] , 

where 

dfXb = da db Q dadddV e- Nia o +b « +M2+v « ) . (4.86) 

Expanding the Fermionic part and keeping only the leading terms in N 2 we 
would find 

u{E) = -^zJd^Al-lal 2 )-"- 2 ^) 2 "- 2 (4.87) 
A [1 - 2iB Q A Q + (iB ) 2 (A 2 a - \a\ 2 ) + 0(1/N)] . 



We use the notations £ = -E/2 + iy/l - E 2 /A so that ££* = 1. 

The ordinary saddle points correspond to \a\ = Vo = and for the main 
saddle point a = ibo = E, so that A = iB = —£* and (iAoBo)^ 1 = (E) 2 . 

At the second saddle point with same bosonic action ao = £,ibo = E* we 
have {AqBq = 1 so that this contribution is 0(1/N). 

Other saddle points have smaller actions £Q. 

The Hessian at the main saddle point is obtained by writing ao — > £ + &o 
and ibg — > £ + ibo (which includes a translation of the integration contour) and 
expanding to second order. The linear part vanishes at the saddle point and we 
find the Hessian approximation: 

v{E) = ^11% f daod6 dada^oe- JV [ (£+ao)2+( - t£+bo)2+|a|2+ ^ 2) ] 

2tt'/ 2 J 

e ~N[\o g [(£+a +E+V a ) 2 -\a\ 2 ]-2\o e (£+ib +E+Vo)] 

£[l-2£ 2 +£ i }+0(l/N) (4.88) 
2^' 

e -N [(a 2 a +b 2 a + \a\ 2 )(l-£ 2 ) + V 2 +2( : ib„-a )V £ 2 ] /jggs 

£[1 -£ 2 } 2 +0(1/N) . (4.90) 

We can compute the Gaussian integral over \a 2 \ as jyfpzgJ) ■ The remaining one 
is (tt/N) 3 / 2 times 1/y/det Q where Q is the 3x3 matrix: 



9 / daodbodadadVo 




(4.91) 
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Since det Q = (1 - £ 2 ) 2 we find: 



v{E) = -$S£ = -v/1-^2/4 



7T 7T 



(4.92) 



which is the Wigner semi-circle law. 

This essentially ends the proof, up to two main effects which have been 
neglected: we have to bound the error to the Hessian approximation to the 
saddle point and to treat the "diagonal" B quartic terms. 

4.3 Bounding the error to Hessian approximation 

First of all, we have to check that the path we have chosen passing through the 
saddle is the correct one. This is true if the real part of the action has minimum 
at the saddle. 

Then we extract the leading contribution from the observable term, (the 
corresponding integral is one by supersymmetry) that gives the semicircle law. 
The remaining error terms are of order 0(1/ N). To prove that we need to divide 
the integration region in the five parameters in four main zones 

- the vicinity of the ordinary saddle point: here the Hessian approximation 
is correct. As we have extracted the leading term already the numerator 



- the vicinity of the second saddle point: also here the Hessian approxima- 
tion is correct, but we have to expand around the second saddle; 

- a large but compact region around the two saddles: here we bound the 
integrand by supnorm; 

- the rest: as we are very far from the saddle we can extract some exponen- 
tial decay. This allows to perform the integral and gives an exponentially 
small factor. 

Extracting the leading term 

Putting together the terms of order N 2 , N and 1 we have 



is small; 



u(E) = 



-Vn 

~A^J 2 



s / <K 



(A 2 -\a\ 2 ) 



A 



[N 2 h +NI 2 + I 3 ] 



(4.93) 



where the measure fj,' b contains also the logarithmic terms: 




(4.94) 
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and 



h 
h 
h 



WTRF) 2 [1 " 2WoAo + {lB ° )2[Al ~ |fl|2)] 

IE 



(iB ) 2 (A 2 -\a\ 2 ) 
1 

~^o) 2 (4MaP) ■ 



l-4A iB + 2 l -^(A 2 -\a\ 2 ) 



(4.95) 



Note that, if we have no observable, the full supersymmetric integral is 1, 
and its approximation to leading order is also 1 ! So we must have 



1 = N< 



47r 5/2 



h ■ 



Therefore 



where 



v(E) = + -9£ + iV 2 i?x + NR 2 + R 3 , 

7T 



(4.96) 
(4.97) 



Ri = 
Rs = 



4tt 7 /2 

-Vn 



4tt 7 /2 

-Vn 



3 / djj! b 



if d^ b 



An 



K - M 2 ) 



h , 



4tt 7 / 2 



S / d v'b , 72— T7TcrT J 3 



(4.98) 



'(^o- N 2 ) 

We need to prove that the three quantities _/V 2 |i?i|, iV|i? 2 | and |i?3| are 0(1/N). 
Translation to the saddle 
We translate to the saddle 



a —> a + £, 6n — > &n - 



(4.99) 



We do not cross any singularity (thanks to the is regulator) hence the integral 
does not change. The measure becomes 



d/i' 6 = da db dadddV e' NS< - aaM - Vo '^ 



(4.100) 



and 



[(£ + a ) 2 + (-i£ + b Q ) 2 + \a\ 2 + V 2 )] 

+ [log[{-£* + a + V f - \a\ 2 } - 2 log{-£* + ib Q + V )] 

[a 2 + b 2 + 2£(a - ib a ) + \a\ 2 + V 2 ] 

+ [log[(-r +a + Vof - \a\ 2 } - 2 log(-£* + ib Q + V )] . (4.101) 
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In order to bound the error terms, we will need a bound on the absolute 
value of the measure. So we have to study the critical points of Re S: 



ReS = a 2 + bl + 2b £ t - 2a a £ r + \a\ 2 + V 2 + - lnF - lnF b 



(4.102) 



where we defined & = y/l - E 2 j^ £ r = E/2, 

F a = x 2 [x 2 + 2(£f-\a\ 2 )] + (£f + \a\ 2 ) 2 



{b + Si) 2 + {V + s r y 



and 



x = a + V + £ r ■ 
The equations for the critical points are 

= d ao RcS = 2 - £ r + ^ [x 2 + £ 2 - \a\ 2 ]^ 



(4.103) 
(4.104) 



= 8 bo ReS = 2(b + Si) 



d a RcS = 2a 



l + U-x 2 +£ 2 + \a\ 2 ) 

f n 



= dy ReS = 2|K,(l--^) + [x 2 +^ 2 -|a| 2 ]|.(4.105) 

We remark that E = 0(X 2 ) = 0(1 /N). Therefore we can take E = in the 
saddle point computation, as it only gives corrections of order 0(l/N). At 
E = we have £ r = and £i = 1. The saddle equations become 

= 9 ao RcS = 2 ja + Jr(a; 2 + 1-M 2 )j , 



= d ba RcS = 2(b Q + 1) 



F b 



d a RcS = 2a 



l + —(-x 2 + l + \a\ 2 ) 

r n 



d Vo ReS = 2 { V 1 - — + — (x z + 1 - la 



F b F a 



(4.106) 



where now x = ao + Vq , 



F a = x 2 [ a; 2 + 2(l-| a | 2 )] + (l + | a | 2 ) 2 
F b = (b + l) 2 + V 2 . 

Note that, from the first and last equation we have 
V 



(l-l)=-J(, 2 + l-|a| 2 )=ac 



(4.107) 



(4.108) 
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Moreover, from the third equation we have a — or (1 — x 2 + \a\ 2 ) — —F a . In 
this last case we have to solve 

x 2 [x 2 + 2(1 - H 2 )] + (1 + H 2 ) 2 = -(1 - x 2 + \a\ 2 ) . (4.109) 

We find there is no real solution for y = x 2 . Therefore we have always a = in 
the following and F a = (1 + x 2 ) 2 . The saddle equations then reduce to 

x 



= a 



2 + l ' 
= (&o + l) 



1 

1 

F b 



- ^( 1 -^) + TT^- ( 411 °) 

Now we distinguish several cases. From the second equation we see that 
b + 1 = or F b = 1. 

First case: F b = 1 

From H4.108fl we have oo = 0. We also have x = Vq, and 

x V 



1 + x 2 1 + Vo 2 



, (4.111) 



which implies Vq = 0. Finally bo is given by 

1 = F b = (b Q + l) 2 which implies b = 0, b = -2 . (4.112) 

Therefore we have two saddle points: do = Vo = a = and bo = 0, —2. 

Second case: bo + 1 = 

In this case we have F b — V 2 , a — V - l/V and x — (2V 2 — 1)/V . 
Inserting the expressions in terms of Vq in the first equation we find 

^1+ V f V »- l) =0 (4.113) 

Vo W a 4 -3V 2 + 1 V j 

which gives the equation 

4z 3 -5z 2 + 3z-l = (4.114) 

where z = V 2 . Studying the first derivative of this expression we see that there 
is only one real solution. This solution is positive with 1/2 < z s < 1. So there 
are two solutions Vo — ±-»/izJ. We have then two new critical points in a = 0, 
bo = -1, Vo = ±V^I and a = (V 2 - l)/V . 

So finally we have identified four critical points: 

51 = a = Vo = a = b = , 

5 2 = a = V = a = 0, b Q = -2 , 

5 3 = a — 0, b Q = — 1, a Q = - " Vo = n/zJ • 



S 4 = a = 0, b = -1, a = -^pS, Vo = -ypT. . (4.115) 
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The real part of the action at the leading order (E = 0) is 

Re E=0 S = a 2 Q + b 2 + 2b + \a\ 2 + V 2 + \^F a - lnF b . (4.116) 

Is is easy to see that in S\ and 52 we have Re_E = oS* = 0. On the other hand in 
S3, S4 the action is (writing z for z s ): 

„ (2z-l)(z-l) 4z 2 -z + l 

Re E =oS=- - ^+log 5 (4.117) 

z z 

which since 1/2 < z < 1 satisfies 

Re E =oS > -1/4 + Iog(3/2) > 0.15 . (4.118) 

Integration regions 

We call X = (ao, 60, Vo, a) so that the integral to study is 

J dX e- NReS W\R(X)\ (4.119) 

where R(X) depends on which error term we are looking at. We introduce the 
norm 

||X|| 2 = X*X = al + bl + H 2 + V 2 . (4.120) 
Now we partition the integration domain in four regions T i: i = 1, ...,4, where 

Ti = \x\ \\X\\ < ' 



]\[l/3+S 



T 2 = {X\\\X-S2\\< 1 



Ni/3+5 

T 3 = {X\ \\X\\ <ifandI^T 1 U T 2 } , 

T 4 = {X\ \\X\\ > K} , (4.121) 

where < 5 « 1 and K » 1 are fixed constants. 
Bound in the region T\ 

In this region we can expand the action to the second order. We get 

iVReS = N[(at + b 2 + \a\ 2 )2£f + V a 2 + 2[b Q 2£ r £ l + a (£i -£?)}V 

+o(\\xf)} 

= NX*HX + 0{N- 3S ) (4.122) 



where 



H = 



( 2£ 2 E 2 -E 2 \ 

2E 2 2£ r £ t 

£ 2 -£ 2 2£ r £ % 1 

V 2£ 2 J 



RcQ 
2£ 2 



(4.123) 
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where Q was defined in i|4.91|) . 

Note that H is positive definite, so the measure is well defined. Now we 
remark that in this region 



h 
h 
h 



(l-£ 2 ) 2 + 0(\\X\\), 
£ 2 (£ 2 -2) + 0(\\X\\) 

-£ 4 + o(\\x\\), 



and 



A 



£ 



.(A 2 -\a\ 2 ) 
Inserting this bounds in Ri,R 2 , -R3 we get 



= 0(11*11). 



(4.124) 
(4.125) 



N 2 \Ri\ Tl < const N 5 / 2 J dXe- NX ' HX 0{\\X\\) =0 
N\R 2 \ Tl < const N^ 2 J dXe' NX ' HX = O (J^j , 
|i? 3 | Tl < constVN J dXe 



-NX'HX _ (~) I 

1 iV 2 



(4.126) 



By parity we can improve the first estimate and obtain in fact rather easily 
N 2 R 1 = 0{\/N). 

Bound in the region T2 

In this region, we expand around the second saddle S2 and we obtain 

NReS T2 = N[(a 2 + (b + 2) 2 + \a\ 2 )2£ 2 + V 2 

+2[(b + 2)2£ r £i + a (£ 2 - £ 2 )]V Q ] + 0(\\X - S 2 \\ 3 )N 
= N(X - S 2 )*H(X - S 2 ) + 0{N- S ) . (4.127) 

The bounds on NR 2 and R3 work as before. For N 2 Ri we have 

h = + O(\\X-S 2 \\) , (4.128) 



A 



£ 



0(\\X-S 2 \\) , 



_(A 2 -\a\ 2 ) 
so that again iV 2 i?i = 0(1/ N). 
Bound in the region T3 

We have to bound f(X) = \ e - NReS(x ) R(X)\. We shall prove 

f(X) = \e- NReS ^R(X)\ < ce- Nl/ " 35 . 



(4.129) 



(4.130) 



The function R(X) is a rational fraction and to bound it we must take care of 
the possible zeroes of the denominator. This denominator is (A 2 , — |a| 2 ) 2 (z£?o) 2 - 
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But (Aq — | a. | 2 ) has no zeroes thanks to the imaginary part iBo in contrast 
vanishes on the submanifold 60 = and Vq = —£r-, but the factor (iBo) 2N 
that we put in S also vanishes. So we are lead to write 

e -"SW R (X) = e -(N-2)s(x)-2S,(x)+^G aR!{x) (4 131) 

where 

G a =(-£* +a + ^o) 2 - |a| 2 (4.132) 
R'(X) — R(X)(iBo) 2 is now a rational fraction with a non-vanishing denomi- 
nator, and 

S 2 (X) = a 2 Q + b 2 + 2£(a Q - ib ) + \a\ 2 + V 2 . (4.133) 
On the compact region K we have 

e -2S 2 (X)+lnG < Kj _ (4134) 

for some constant K\, and 

\R'{x)\<cN 5/2 (l + \\X\\) d (4.135) 

for some constants c and d. 

So it is sufficient to prove that in that region 

(N - 2)ReS T3 > c(N - 2)^^ - c'N 1 ^ 26 (4.136) 

for some small positive constant c and d . Indeed we have to check where KeS 
can reach its minimum. This can only be on a critical point or on the boundary. 
But the only critical points in T3 are S3 and S4, where we have the much stronger 
bound NKeS > 0.15AT. On the boundary with T\ and T2 the previous Hessians 
approximations are still valid and give the bound 1)4.1361) . Finally on the outer 
boundary the bound is even better (see below). 

Bound in the region X4 

When we are far enough from all saddle points we can apply a rough bound 
to the logarithms in the action 

e A'[ilnF a -(A'-2)lnF b ] < e +K'N(lo S \\X\\) (4.137) 

for K' big. Moreover 

e +J*r'JV(iog||x||) < e +N\\xf/A (4.138) 

taking K (the parameter defining region T4) large enough, since recall that in 
Ti \\X\\ > K. 

Finally in T4 it is easy to check that for large enough K 
(N - 2){a 2 + b 2 c + 2b £ t - 2a £ r + \a\ 2 + V 2 ) > (l/2)N\\X\\ 2 (4.139) 
so that remembering l|4.102|l and the bound 14.135|) on R' we can bound 

/ dXe- NRcS ^\R(X)\ < [ e - NVertx ^/ 4 N^ 2 K 2 (l + \\X\\) d < e~ cN . 

Jt 4 J\\x\\>k 

(4.140) 

This completes the proof. 
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4.4 Bounding the correction due to B 

This is explained in a rather detailed way in the last section 6, of PP, which for 
completeness we roughly reproduce here. 

The quartic term B should be small as 1 /N compared to the other ones since 
it contains only one sum over sectors a, not two independent sums over a and (3. 
However it has the wrong sign, so although statistically small, it is dangerous 
at large fields. The solution is to treat it by some kind of small/large field 
expansion. A Taylor expansion with integral remainder is written successively 
for each of the N sectors appearing in the sum for B: 



B = Y J B a . (4.141) 

a 

To first order the expansion gives 

l 

e- x2B « = 1 - / \ 2 B a e- x2B «e +tx2B « dt. (4.142) 



This Taylor expansion either suppresses B a from the exponential of the action 
or generates a remainder term 



R a = - I \ z B a e- x B « e +tx B « dt. (4.143) 



The rough line of argument is as follows: either the set P of sectors where B a 
is not suppressed is small or it is large. Let p — \P\. If p is large, we have 
gained such a small probabilistic factor that we can treat this case by a rough 
translation of the V contour, like for large field cubes of section (3.3) (regime 
\p 2 — E\ » A 2 ). No saddle point analysis is necessary. 

If p is small, we have to perform the saddle point analysis, but we can restrict 
is to the sum of sectors outside P where it is essentially the same analysis than 
previously but with only N' = N — p sectors. We have to treat also the coupling 
between the sectors in and outside P, and this is a delicate point, since the 
ultimate reason for which these couplings are small is the supersymmetry of the 
underlying field theory. 

Now what " small" or " large" means for pi It means p smaller or larger than 
some function p(N). Since we expect a small factor of order 1/iV for each sector 
in P (through the lacking sum over /3), and a rough imaginary translation on 
V costs K N , we need K N e- p ^ losN « 1. 

Therefore we choose for p(N) the integer part of N / \f\og N . The expansion 
is stopped at order p{N) = p. This means that we write: 

e- x2B = l + ^Pc[l, N]0 < \P\ < p(N)R P , (4.144) 
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where 



r p = n r « if i p i < ?w. 

R P =l[R a J] e- A2e ° if|P|=p(JV). (4.145) 

a £ P a > max P 

Let Q be the complement of P in [1, ...,JV]. The term 1 was treated in the 
previous sections. The remainders terms Rp must be shown to be 0(1 /N) as 
N -> oo. Let i? p(Ar) = E p C [1,...,7V]|P| =p(N)R P be considered first. For 

this term we said that it is not necessary to perform any saddle point analysis. 
It is sufficient to return to the treatment of subsection 3.3, so wc undo the V 
integration: all the e~ A Ba -terms are recombined with the A-teim to reproduce 
the initial functional integral 1)4.72)1 . However we cannot exactly reproduce 
the initial functional integral over V, since some e _A 8 ° -factors are missing or 
appear with reduced weights in (|4. 145(1 . This means that there remains quartic 
correction terms e tx Ba or e A Ba . The important remark is that the bosonic 
part of these terms has now the right sign ! Therefore they can be represented 
as a well defined functional integral over a new auxiliary field W a . For instance 

exp ( - 2\ 2 tS+S- a Sl a S a ) = cxp ( - 2X 2 t\S a S- a \ 2 ) 

= J dW a dW a e - |H/ ° |2 exp (i\V2tW a S a S^ a + cc). 

With slightly condensed notations, this leads to 



Rp(N) 



E / S+S (U J a2 " 



PC[1,...,JV] , |P|=p(Af) 1 ' aEP'j 

= x cxp (m + (E + XV + \V2tW)^) dn(V, W). (4.146) 

Then a complex translation Vo 1 — *• Vb ± iX^ 1 is performed, with the same sign 
as the imaginary part of E in order to avoid crossing of singularities. In other 
words 

e- v °F(V )dVo = { e- v o e - 2lX ~ lv °e x ~ 2 F(V +iX- 1 )dV . (4.147) 



The functional integral can now be bounded by its absolute values everywhere, 
namely the following contributions are bounded 

- by 2 N , for the sum over P, that is the total number of subsets of [1, N]; 

- by 1, for the integrals such as J Q dt; 

- by 1, for the oscillating imaginary integrals; 
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- by Gram's inequality for fermions or the Schwarz inequality for Bosons, 
for the remainders terms B a ; 

- by 1, for every propagator since the imaginary translation in Vq has created 
an imaginary part proportional to the identity in the denominator of the 
Green's function. 

This means that each B a term gives rise as expected to a small factor A 2 = 1/2N 
for each sector a in P, hence altogether we have a small factor 1 /N p ^ . The two 
source terms are bounded by 1. The normalization determinants arc then easily 
bounded by c N , even without using the supersymmetry cancellations, since the 
operators considered are bounded in a finite 2N dimensional space thanks to 
the imaginary part of E which is no longer infinitesimal. Combining all factors 
leads to 

R P( N) < c^e-^, (4.148) 

showing that this correction term is indeed small. 

It remains to treat the terms with 1 < |P| < p(N). To bound these terms 
a mean-field analysis will be performed like in subsection (3.2), introducing the 
a and b fields, but we said it should apply only to the sectors of the theory in 
the complement Q of P. The functional integral to be bounded for a single 
term is <|4.145[1 . Now, the e~ x Ba -terms are recombined only for a G P with 
the A term to reproduce the initial functional integrals over the V fields and 
the terms e tx BoL , also for a £ P, are again given by defined integrals over new 
auxiliary fields W a . Finally, the quartic terms, with sector sums reduced to 
Q, are treated exactly as in the previous section, hence mean-fields a, ao, &o & r e 
correspondingly introduced. This leads to a representation 

R P = I dV Pt p dV P , Q d 2 a dV da Q db J| R a e CQ , (4.149) 

J a£P 

where Q is the complement of P in [1, N]. In addition, Vp.p is the part of the 
matrix V corresponding to rows and columns in P, including the new fields of 
the W type. Vp,q correspond to one entry in P and the other in Q, and the mean 
field computation is now restricted to Q. The integral over superfields gives rise 
again to a superdeterminant and an additional correction, whose bosonic part 
(leaving the Fermionic part to the reader) is of the type 

exp[Trlog(l + CV)], (4.150) 

where C is the bosonic part of the matrix for the N — p analogous problem, 
as in the previous subsection (see eq. l[4.78Jl ). and V = Vpp + Vp,q is the 
perturbation. This correction term is bounded by 

| exp [Tr log( 1 + K)] \ < cxp(TrK + K* +KK*). (4.151) 

Evaluating C — R^ 1 at the saddle point costs a factor exp (Np/2N) = e 2p 
at most. Each term B a naively gives a factor A 2 = 1/(2N) when evaluated, 
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but this simply compensates the sum over a when p is nonzero but small, so 
we have to gain an additional 1/N factor. Adding a few expansion steps gives 
such a small additional factor 1/N (already for the first non zero value p = 1). 
This can be seen by integrating by parts the superfields in the vertex B a which 
has been taken down the exponential. By supersymmetry, the vacuum graph 
corresponding to a contraction of the four fields at the vertex vanishes. This is 
absolutely necessary since this graph by simple scaling is proportional 1 /N and 
cannot have any additional 1/N factor. Its vanishing can be checked by hand: 

- the self-contractions of the bosonic piece 2S^ S- a S^ a S a give a factor +2; 

- the self-contractions of the boson-fermion piece 



give a factor -2 (since there is one Fermionic loop giving the minus sign 
and one sum over a giving a factor 2 only after the contractions); 

- the selfcontractions of the term (1/2) ^^i^ta^cra) 2 are clearly supersym- 
metric and also add up to 0. 

Consequently, at least one field of the vertex B a has to contract to the exponen- 
tial. Performing two contractions in turn, gives always at least a factor p/N 2 at 
the end instead of the naive 1/N factor. Indeed the worst case corresponds to 
the non trivial contraction term being of the type Vpq . This generates a new 
factor 1/N but a new sum over (3 £ Q, which costs N — p ~ N so nothing is 
gained yet. But contracting this (3 field again either generates a diagonal term, 
hence a new 1/N factor, with no new sum, or returns to a Vq p term. This last 
situation generates a new 1/N factor and a new sum over sectors 7 but this 
time this new sum is restricted to P, so it costs only a factor p instead of AH 
Hence at worst, after these two contraction steps, a total factor p/N 2 instead 
of the naive factor 1/N is associated to each vertex, as announced. Now the 
sum over P costs a total factor N\/p\(N — p)\, hence is bounded by c p [N/p] p . 
Combining all factors, the sum of contributions of such terms with p 7^ is 
bounded by J2p=ii c /N} p < d/N. It is therefore at least as small as 1/N. 

This completes our sketch of the content of The last section in this 
review is devoted to generalizations that better mimic the symmetries of the 
2-dimcnsional matrix model near degeneracy of the rhombus, using anisotropic 
sectors. 

5 Improved Flip Matrix Models 

As we already said in Section 3.2, the isotropic sector model is not optimal for 
momentum conservation analysis. 

Actually we know that the potential V(j), has to be associated to pairs of 
vectors p, q with \p\ = 1 + 0(M _J ), \q\ = 1 + 0(M _J ) (using the same notation 
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as in previous sections) and p — q = 7. By simple geometric arguments we can 
see that the angle a between p and q is then sin(cr/2) ~ M/2 and p can move 
of an angle 8a = 0(M~ 3 /I7Q. 

- If we have isotropic sectors (of width M~ J ), then for each |-y | there is a 
set 5(7) of 0(1/|7|) pairs of sectors a,f3 (with angle a between them) 
associated to the random variable ^(7). So when |-y | < 1 the matrix ele- 
ments V a p for (a, f3) € S(j) are no longer independent random variables. 
In the limit case, I7I = 0{M~ 3 ) all the 0(M 3 ) sectors are associated to 
the same V(j). This is the diagonal term Vo that appears in the Isotropic 
flip model. 

- If we still want to work with matrix elements which are all independent 
random variables we have to introduce anisotropic sectors of width M~ J / 2 . 
In this case the matrix elements are all independent for 7 > M~ 3 / 2 . The 
case 7 < M~ J / 2 corresponds to the diagonal terms V aa . Here we do 
expect to have no longer independence, but the problem arises only on 
the diagonal terms, just as in the Isotropic flip model. The drawback 
is that now, for each sector a we have also to introduce a dual space 
decomposition R a of the spatial cube A of side M 3 in rectangles of side 
M J in the direction parallel to a and M 3 ! 2 in the orthogonal direction. 
The random matrix is then indexed by pairs (a, t) where a is a sector and 
t a dual rectangle. 

5.1 The Anisotropic (2d) Flip Matrix Model 

As in the previous sections we consider the Anderson model in a spatial cube A 
of side M 3 , and we introduce a cut-off in momentum space that selects only an 
annulus of width M~ J around the Fermi surface. We fix then M 3 = A~ 2 = 2N, 
E = 1 and E = 0{\ 2 ) = 0{N- r ). 

We then cut the annulus in anisotropic sectors of width M~ J / 2 , so that a 
sector now has width M~ 3 in the radial direction and M~ J / 2 in the tangential 
direction. Therefore we have M J / 2 anisotropic sectors instead than M 3 isotropic 
ones. For each sector a then we introduce a partition R a of the cube A in 
rectangles of side M J in the direction parallel to a and M J ' 2 in the orthogonal 
direction. 

In this way we obtain a finite dimensional space of dimension 2N = M° with 
a new index for its orthogonal basis, namely the pairs (a, r, a) where a runs 
over the first V^V anisotropic directions, r € R a runs over the y/W rectangles 
elements of R a and as previously a takes two values ±1 corresponding to the 
two opposite sectors on the circle. The dimension of the Hilbcrt space for the 
matrix model we want to construct is then the same as in the isotropic case. 

We list now the columns and rows of our random matrix as before, with 
sectors from 1 to V^V, then from —\/~N to -1, and for each sector a we list the 
new rectangular variables r in direct order according to an arbitrary ordering 
of R a - Note that the rectangles for a and —a are then the same modulo a 
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rotation of ir. The effect of the rotation is that the ordering of the rectangle 
is inverted (see FigEJ). This is actually good as it ensures the symmetry of the 
matrix around the anti-diagonal. 

Finally, as in the case of the isotropic sectors, we discretize the potential 
and build a random matrix model, that will be called Anisotropic Flip Matrix 
Model (AFMM). 

5.1.1 Construction of the matrix model 

To construct the matrix we introduce a discretization of the phase space for V, 
in such a way to ensure that: 

- the off-diagonal matrix elements are independent random variables with 
covariance one. 

- the dependence among the diagonal terms is explicitely extracted. This 
means each V aa is a sum of independent random variables V aa — Vq + 
Vi + ... + Vj/2 where Vq appears in all diagonal terms and Vj/ 2 only in one. 

Note that as V is restricted to act on the shell \p 2 — 1| = 0(M~i), the 
momentum 7 for V is restricted to I7I < 2 + 0(M~°). In position space V is 
restricted to the cube A. 

To construct the off-diagonal matrix elements we divide the momentum re- 
gion 7 in three parts: 

• a forward-momentum region corresponding to M~ • ? / 2 < I7I < 1 (sector 
pairs at an angle M _J / 2 < a < 7r/3), 

• a backward momentum region corresponding to 1 < |7| < 2 + 0{M~ 3 ) 
(sector pairs at an angle 7r/3 < a < ir + 0(M~ J )), 

• a 0-momentum region corresponding to \p\ < M~ 3 / 2 (a < M _J / 2 ); this 
last region gives the diagonal terms of the matrix 

In these regions the discretization of the phase space for V must be different. 
To see that let us consider a pair of sectors a, [3 with two vectors inside p£ a, q € 
(3 and 7 = a. — (3. As p and q move inside their sectors the variation for 7 is 
given by 

(57,1 = NT 3 ' 2 cosQ + M- j/2 ) Sj ± = Af^' /2 sin(| + M~ j/2 ) . (5.152) 
as M~H 2 is the sector size. Now we distinguish four limit situations (see Fig l5IGH 

1. When |7| ~ 1 (a ~ tt/3) then Sj\\ ~ 5~fj_ = 0(A/" 3 / 2 ) so we should cut 
momentum space in cubes of side M~ 3 1 2 . This corresponds to 72 in Fig|SJ 
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Figure 5: Variation of 7 as p and q move inside their sectors 



2. When | 7 | ~ M"* (a = 0(M^'/ 2 )) we are almost on the diagonal part 
of the matrix (71 in Fig|SJ. The value of 7 is fixed with a precision much 
higher than M~H 2 in the perpendicular direction, that is <57|| = 0(M~ J / 2 ) 
while 5t j_ = 0(M~ J ). In this region we must then use rectangles of width 
M~3' 2 in the radial direction (that of 7y) and Af~ J in the orthogonal direction. 

3. When | 7 | = 2 — Af~ J (c ~ it) the situation is reversed. The vector 7 is 
fixed with high precision in the parallel direction. For <57y we have 

*7|| = Ar j/2 y/l - (|7|/2) 2 = M- j/2 y/M-i - M-V/4: = 0(M~ j ) (5.153) 

while ^7_l = 0(M^' 2 ). So this time we must use rectangles of width M - - 7 in 
the radial direction (that of 711) and M~M 2 in the orthogonal direction (see 73 
in Fig El . 

4. Finally when | 7 | < M~ j / 2 (a < M~ j ' 2 ) we are on the diagonal terms of the 
matrix. Here too we have two extreme cases. For | 7 | ~ M~^ 2 the momentum 
space should be partitioned as above, and there is still only one sector pair (so 
only one matrix element V aa associated to the corresponding random variable. 
But when | 7 | < M~ 3 , then we have only one cube of side M - -? around the value 
7 = and the corresponding random variable appears in every diagonal term. 

In order to interpolate between these extreme situations we have to introduce 
shells where I7I is approximately constant. In each shell then we introduce cubes 
or rectangles depending on the value of | 7 |. The details are shown below. 

Forward- momentum. This region corresponds to the annulus between I71 1 
and I72I in FigEJand is divided into shells according to an index k, k = 0, ...,j/2— 
1 so that M~ k - X < \p\ < M~ k . 

Each forward-shell is now further divided into rectangular momentum cells 
of length 0(M _J ' 2 ) in the direction parallel to the center of the cell, and of 
width M~ k ~H 2 in the direction perpendicular to the center of the cell. 
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Finally, for each such cell, the cube A is divided into corresponding dual 
tubes t of length M k+ ^ 2 in the direction perpendicular to the center of the 
cell, and of width M^ 2 in the direction parallel to the center of the cell. Our 
normalization is that each random variable for a pair (7, r), which corresponds 
to a volume 1 in phase space has covariance 1. These variables are independent 
Gaussian variables. 

Backward momentum region. The backward region (region between I72I 
and 1 73 1 in Fig[H|) is also divided into shells labeled by a similar index k = 
0, j/2, but the backward fc-shell corresponds now to momenta 2 — M~ 2k < 
\p\ < 2-M- 2k - 2 for k < j/2. The last shell is simply 2-M~i < \p\ < 2 + M~^. 
Each shell is now further partitioned into 0(M^~ k ) cells 7 of width M -J / 2 in 
the direction perpendicular to the center of the cell, and length M _J ' 2_fe in 
the direction parallel to the center of the cell. The cube A is again divided into 
dual tubes r of width M^ 2 in the direction perpendicular to the center of 7 and 
length M^ 2+k in the parallel direction. Each corresponding random variable 
Vy.r corresponds again to a different discretized independent random variable 
for the potential V with covariance 1 since it corresponds to a volume 1 in phase 
space. 

0-momentum region. The last shell is simply \p\ < M~^ 2 . The treatment 
of this region is subtle because the anisotropic sectors are not a basis very 
well adapted to this zero momentum region. We remember that the coefficients 
cannot be totally independent because in the isotropic basis (treated previously) 
there is a diagonal term proportional to a single variable Void where Id is the 
identity matrix, and this term must be present in any basis. They are not 
totally dependent either, so the best description consist in further splitting 
the phase-space corresponding to the 0-momentum region into unit volumes in 
phase-space. 

Therefore we split the 0-momentum region into shells Sh Q ^, in which the 
momentum size is between M~i +k ~ 1 and M~i +k (for k' = the lower bound 
being simply 0). Each such shell is divided into 0(M k ) cells 7' of length 
0(M _:,+fe ) in the direction parallel to the center of the cell, and M~i in the 
perpendicular direction. The central shell k' = corresponds to a single cell 
(the "true" zero momentum region). 

The size of the cells ensures that, given a shell k! and a matrix element V aa 
only one cell in the shell contributes to it. On the other hand the same cell may 
contribute to different diagonal terms. In the limit k' = the same cell appear 
in all diagonal terms. 

For each cell we divide the cube into dual tubes r', with length NP in the 
direction perpendicular to the center of the cell, and width IvP~ k in the parallel 
direction. Remark that the total number of pairs (7', r') which pave the phase 
space of this 0-momentum region remains 0{M : >) as it should be, and that to 
each such pair is associated a random variable of covariance 1. 
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Figure 6: Different decomposition depending on | | : for |^| = I73I ~ 2 we have 
to use rectangles of side M ~ J in one direction and M _J / 2 in the other direction, 
for I7I = I72I ~ 1 we have to use cubes of side M~ J ' 2 and for I7I = I71I ~ M~ J ' 2 
we have to use again rectangles, but exchanging the sides 

The resulting random matrix. Now to model the rhombus rule we have 
to understand the random variable V coupling two anisotropic pairs (a, t) and 
(f3, t'). A momentum cell 7 of V defines uniquely four possible pairs a, j3 (up to 
nearest neighbors, neglected). Indeed we still have Hermiticity and spin flip: 

V(a, 13) = V(a, 0) = V(-0, -a) = V(-a, -(3) . (5.154) 

Accordingly there are about 0{M k ) pairs a, (3 which have their sum in the 
forward or backward shell of index k. We consider then two tubes t, t' and 
remark that their intersection is either empty or roug hly of width M J/2 and 
length A / P/ 2+k . Furthermore as soon as k increases, the direction of a and (3 
begins to collapse. This applies both to forward and backward scattering. 

When they intersect, their intersection is a tube of length M k+ ^ 2 and width 
M J / 2 and there is therefore a single independent variable V ltT of covariance 1 
both in the forward scattering case and in the backward scattering case corre- 
sponding to this intersection (see Figur^7|). 

When the two tubes miss (see Figure |SJ, we put a zero in the corresponding 
random matrix. 

This describes the random matrix if we complete by Hermiticity and sym- 
metry around antidiagonal (spin flip), except for the diagonal blocks (a, t; a, t') 
which correspond to the 0-momentum region for V . The coefficients of these 
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blocks are except on the pure diagonal t = t', since two distinct parallel tubes 
t and t' are disjoint. (Similarly the antidiagonal blocks corresponding to ex- 
treme backward scattering a = —(3 are purely antidiagonal (0 except on the 
anti-diagonal). But the antidiagonal is made of independent random variables 
of covariance 1, while the diagonal is more complicated). As we have seen above, 
a diagonal term is represented as a sum of independent random variables each 
corresponding to a different shell k' . 

Let (a,t) be a given anisotropic sector and tube, and k' a shell index. We 
call -f'(a,k') the cell of the shell of index k' perpendicular to the center of a, 
and we call T r (t,k') the corresponding associated tube containing t. The pair 
(7' (a, k'), r'(t, k')) corresponds to a single random variable of scale k' associated 
to (a,t). 

The diagonal coefficient of our random matrix is then simply 

M(a,t);(at) = ^ V~j' {a,k>) ,t' (t,k') . (5.155) 
fe'=0 

Remark that conversely for a given {'y',T / ) in shell k' there are roughly 
Mj/2-k tubes t in r', and roughly M^ 2 ~ k cells a giving the same 7', hence 
M J_2fc pairs (a, t) giving the same {'y' ,t'). This means that a given random 
variable Vy, r ' appears at M^~ 2k different diagonal places in the matrix. For 
k' = the single random variable, also called Vo appears everywhere. At the 
other end, for kl — j/2 the variables Vy/,T' appear essentially at a single place. 

5.2 The hierarchical approximation 

This approximation consists in somewhat simplifying the intersection rule for t 
and t' and the angular condition on a and (3. For (a ^ j3) instead of considering 
that the pair is of class k if their sum lies in the fc-th shell, we divide the set of 
anisotropic sectors into M k packets for each value of k, k — 0, ...1, ■■■j/2 (each 
packet containing M^ 2 ~ k sectors) and we replace the rule on the sum 7 to be 
in shell k by the simpler "hierarchical" condition that a and (3 belong to the 
same packet of scale k but not to the same of scale k + 1. 

In the same way we simplify the intersection rule for the tubes. We divide 
for each sector a the set of tubes T a into M k packets, according to translation 

(k) 

in the direction perpendicular to a and label these packets with an index r al — 
1, M k . Each packet corresponds itself to a larger tube of width M J_fc , hence 
contains M^ 2 ~ k anisotropic tubes of width M^ 2 . Let us consider now two 
tubes t, t' for a pair (tr a, erg/3) of category k, where a a = ±1 says if a is in the 
upper half or lower half of the circle. Then we slightly simplify the intersection 
rule by saying that, if a a = ap, then t, t' intersect if their packet index of degree 
k is the same: r4 — and miss completely otherwise. On the other hand, 
when a a — — ap, for instance a a = l,ap = — 1, then the packets in Rp have 
performed a rotation by an angle 7r, and the ordering is reversed (see FigEJ. In 
this case we say that t, t intersect if their packet index satisfy r a = M —r\ . 
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Finally we take M = 2, which also simplifies the computation. In the limit 
N — > oo the density of states of the real problem should be similar to the limit 
of this hierarchical limit. 

The matrix can then be represented as a set of %/iV x y/N blocks. Each 
block corresponds to a fixed sector pair, and the variables in the block give the 
couplings between the tubes. In the blocks corresponding to k — 1 (both in the 
backward and forward region) all sectors intersect, so all elements are indepen- 
dent random variables. For k > 1 the block has a block diagonal structure, with 
2i/2-fe x 2i/2-' c blocks along the diagonal or the antidiagonal, depending if we 
are considering the forward or the backward region. The elements inside the 
blocks are independent random variables, the elements outside are zero. In the 
limit when k — j /2 the blocks have non zero elements only on the diagonal or 
antidiagonal. In the first case, the elements are no longer independent variables. 

Note that, from the definition of forward and backward shells, the two k = 
shells cover angles tt/6 < a < 2/3ir approximately. In this region k = then all 
tube packets intersect and the corresponding block of the matrix has no zeros. 
Therefore, for each line in the matrix, at least half of the blocks have no zeros. 
For k > (both in the backward and forward region) we are considering only 
small angles. As a result the matrix looks like FigllOl 

We still have to say what are the diagonal terms V a .t- a ,t- We may consider 
only the first \/N sectors as the others are obtained by flip symmetry. The 
diagonal T> can then be written as a sum of diagonal matrices 

J/2 

V= s £ j V k , . (5.156) 

k'=0 

For k' ^ k" the corresponding random variables are independent. For a fixed k' 
T>k' is again a diagonal matrix, made of 2 k diagonal blocks B(k') r r = 1, ...2 k . 
Each block B(k') r actually selects a subset of 2^ 2 ~ k sectors that are associated 
to the same set of random variables. For kl — there is only one block as 
2?o = Void. In the opposite case k! — j/2 there are 2 k as all diagonal elements 
are independent random variables. As an example, for k' = 1 we have only two 
blocks 

* - ( T ) ■ 

Then a single block B(k') r is built from 2 J / 2 ~ fe identical sub-blocks V(k', r), 
one for each sector. This sub-block V(k',r) actually identifies the set of inde- 
pendent random variables associated to this set of sectors. For k! = the single 
block -B(O) has 2^ 2 sub-blocks, that is all sectors are associated to the same set 
of variables. 

Finally, we have to write the y/~N x V~N block V(k', r) which corresponds 
to the tube indices. This block contains 2 k independent random variables. In 
order to distribute them on the diagonal we introduce a partition of the set 

{1, \ZJV} in 2 k ' sets of 2^ 2 - k ' indices, defined by / = Llf^Jj and 

h={i\ 2 j ^- k \l -l)<i< 2^ 2 - k 'l} . (5.158) 
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Figure 10: Form of the random matrix after the hierarchical approximation; the 
F represent the areas where no zeros appear, D and A are the areas where only 
blocks on the diagonal or antidiagonal are non zero. The size of the block is 
very big as we are near to the border with the F region and reduce to one as we 
go towards the diagonal and antidiagonal 



t/6 71/2 




Figure 11: Form of the random matrix if we select only the first /2 sectors 
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Now to each I = 1, 2 k we associate an independent random variable V(k' , r)j. 
With this notations 

V ii = V{k',r w ) l , h ifieJj. (5.159) 

Again we see that in the case fc' — we have only one set I so that there is only 
one independent random variable. 



5.2.1 Toy model 

To understand how the analysis is going to change with respect to the isotropic 
model we first consider a simplified model, where we restrict to the sectors 
1, ...\f~N /2, that is to angles < a < n/2. The matrix then looks like FiJTTI 

We can construct now this matrix in a hierarchical way. 

The first hierarchical step is fc = 1, with two packets of sectors, and two 
packets of tubes. If we list the indices of the first packet first and second packet 
second, and then within each sector packet we list the tube indices with first 
packet of tubes first, this cuts the matrix into 16 blocks. Then we put random 
independent coefficient (with the Hermiticity condition) everywhere except in 
the two upper blocks Bn and #34 which we fill with zeros. (We complete 
by Hermiticity so that lower blocks #21 and #43 are also zero). These zeros 
correspond to the first group of missing tubes. In this way we get an N x N 
called matrix Mi (N) 



M X {N) = 



( #11 





B21 


#13 
#23 


#14 
#24 






(BO 
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B* 




#23 
#24 


#33 






#44 
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B 


B 
B 



where B^ = Ba for alii = 1, .., 4 and each block Bij is a N/4 x N/4 matrix. 

Then for k = 2, the 4 previous diagonal blocks of Mi(iV) are themselves 
replaced by Mi(N/4). This defines M2(N), which is a sixteen by sixteen matrix. 
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We iterate. To get Mk(N) we replace each of the 4 fc 1 diagonal blocks of 
Mk-i(N) by Mi(N/A k ~ 1 ). We stop at k = j/2. Then the diagonals blocks of 
size 1 have been reached. 

Finally in the last step we change the diagonal coefficients of Mj/% (N) . We 
replace each of them by the appropriate sum of independent variables according 
to the analysis of the zero momentum region. Again this is done in a hierarchical 
way. At k = we put a Vq at every diagonal point. At k = 1 we split the 
diagonal in two parts. To all terms Vu with i < N/2 we add the variable Vy, to 
all others Vf. To iterate that we assign an ordering r-fc(i). For k = 1 we have 
fk{i) — 1 for all i < N/2 and rk{i) = 2 otherwise. At k = 2 we split in four sets 
and iterate. 

The diagonal coefficient V^i of the final matrix is therefore 

V u = V + ... + Vl k{€) + V]]^ (5.161) 

where at each scale k there is 4 fe ~ 1 variables VI, each serving for N/4 k ^ 1 coef- 
ficients, and rk(i) — 1, 4 fc_1 is simply the label for the variable of the cell in 
the shell k associated to i. 

The rest of the matrix remains unmodified, with independent variables. 

So the matrix is made of 

- independent random variables at each non-diagonal place, except where 
the zeros of Mj/ 2 (N), lie, which correspond to tubes missing each other, 

- a more complicated sum of random variables at each diagonal place. Re- 
mark in particular that this sum always contains the Vq variables, common to 
all the diagonal. 

Again the density of states for this matrix should follow some semi-ellipse 
Wigner's law when N — > oo. We have to reproduce the computation of the 
previous section but with additional fields which complicate the picture. Each 
set of zeros should indeed be represented by a new pair of complex fields which 
suppress the off diagonal blocks with zeros in the quartic form after integrations 
of V. 

When searching for a main saddle point, an iterative solution seems doable 
by plugging the smaller blocks solutions into large means fields, until the last 
one for which we should get v{E) = Ck \/ E'i — E 2 where Ek lies between y/2 
and V3 GUI- 

5.3 A Hierarchical Isotropic matrix with correct Momen- 
tum conservation rule 

An alternative is to stick to the initial basis of isotropic sectors but to take 
seriously into account the necessary modification of the rhombus rules for de- 
generate rhombuses. 

The approximation in which each random variable corresponds roughly to a 
single pair of momenta, up to Hermiticity and flip is indeed not correct. When 
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the momentum of V is of size M~ k (forward) or 2 — M~ 2k (backward), there 
are roughly M k pairs that have this transfer momentum. 

Taking again M = 2 the hierarchical approximation to this correct conserva- 
tion rule consists in identifying blocks of 2 fe by 2 fc coefficients to a single random 
variable as approaching the diagonal (or anti-diagonal in the version with the 
flip symmetry). 

In this point of view the hierarchy is indexed by k from to j (not j/2). 
This model should again give the same limit as the previous one when j — > oo. 

6 The isotropic 3D Model 

We would like to conclude with a remark concerning the three dimensional 
model. The U(l) symmetry which generalizes the flip symmetry in that case 
has orbits which are quite complicated to describe because there is no simple 
labeling of sectors on the sphere which respect the rotation invariance of the 
sphere. For instance labeling according to spherical coordinates introduces a 
preferred polar axis. So here again it may be instructive to mimic this U(l) 
symmetry by treating first some hierarchical approximations in which the U(l) 
symmetry is replaced by a discrete Z p symmetry, with p = 2 k . For instance 
such symmetries could act in a simple way on coefficients, by folding k times 
the upper triangle around its bisecting line, and identifying the 2 k coefficients 
at the same position in all the layers. Then letting p ~ ^/N and N — > oo, it 
should be possible to prove again through a supersymmetric analysis that the 
density of states of such a model tends to Wigner's law. 

Such a result would certainly increase our confidence that the density of 
states for the true three dimensional Anderson model falls into the same uni- 
versality class than the GUE and obeys Wigner's semi-circle law. 

References 

[1] J. Bcllissard, Jacques Magncn, V. Rivasseau, Markov Processes and Re- 
lated Fields, 9, 1-30 (2003) 

[2] D. Voiculescu, Limit Laws for Random Matrices and Free Products, Invent. 
Math. 104, 201 (1991) 

[3] M. Disertori, H. Pinson and T. Spencer, Density of states for random band 
matrix, Commun. Math. Phys, 232, 83 (2002). 

[4] M. L. Mehta, Random Matrices, Academic Press, 1991 

[5] L. Pastur and A. Figotin, Spectra of Random and Almost-Periodic Oper- 
ators, Springer- Verlag 1982 

[6] A. Mirlin, Statistics of energy levels and eigenfunctions in disordered and 
chaotic systems: supersymmetry approach. Phys. Rep. 326, 259-382 (200) 



50 



[7] M. Zirnbauer, Journ. Math. Phys. 37, 4986 (1996) 
[8] G. Poirot, These, Ecole Polytechnique, 1998 

[9] B.I. Shklovskii and A.L. Efros, Electronic Properties of Doped Semiconduc- 
tors, Springer (1984) 

[10] P.W. Anderson, Absence of diffusion in certain random lattices Phys. Rev. 
109, 1492-1505, (1958) 

[11] E. Abrahams, P.W. Anderson, D.C. Licciardello and T.V. Ramakrishnan, 
Scaling theory of localization: absence of quantum diffusion in two dimen- 
sions. Phys. Rev. Lett. 42, 673-676, (1979) 

[12] I. Gold'shcid, S. Molchanov and L. Pastur, A random homogeneous 
Schrodingcr operator has a pure point spectrum, Fund. Anal. Appl. 11, 
1-10 (1977) 

[13] H. Kunz and B. Souillard, Sur le spectre des operateurs aux differences 
finies aleatoires. Commun. Math. Phys. 78, 201-246 (1980) 

[14] J. Frohlich, F. Martinelli, E. Scoppola and T. Spencer, Constructive proof 
of localization in the Anderson tight binding model. Commun. Math. Phys. 
101, 21-46 (1985) 

[15] M. Aizenman and S. Molchanov, Localization at large disorder and at ex- 
treme energies: an elementary derivation, Commun. Math. Phys. 157(2), 
245-278 (1993) 

[16] G. Bcrgmann, Weak localization in thin films, a time-of-flight experiment 
with conduction. Phys. Rep. 107, 1-58 (1984) 

[17] B. Kramer and A. MacKinnon, Localization: theory and experiments. Rep. 
Prog. Phys. 56, 1469-1564 (1993) 

[18] K.B. Efetov, Supersymmetry and theory of disordered metals. Adv. Phys. 
32, 53-127. (1983) 

[19] I.Kh. Zharekeshev, M. Batsch and B. Kramer, Crossover of level statistics 
between strong and weak localization in two dimensions. Europhys. Lett. 
34, 587-592, (1996) 

[20] B. Altshuler and B.I. Shklovskii, Repulsion of energy levels and conductiv- 
ity of small metal samples. Sov. Phys. JETP 64, 127-135 (1986) 

[21] K.B. Efetov, Supersymmetry in Disorder and Chaos. Cambridge University 
Press, 1997 

[22] Jean-Michel Combes, Peter D. Hislop, Frederic Klopp and Georgi Raikov, 
Global continuity of the integrated density of states for random Landau 
Hamiltonians, Proceedings of the Marseille conference in honor of L. Pastur, 
2002 



51 



[23] G. Benfatto and G. Gallavotti, Perturbation theory of the Fermi surface in 
a quantum liquid, a general quasi-particle formalism and one dimensional 
systems. J. Stat. Phys. 59, 541-664 (1990) 

[24] J. Fcldman and E. Trubowitz, Perturbation theory for many fermions sys- 
tems. Helv. Phys. Acta 63, 156-260 (1990) 

[25] J. Feldman, J. Magncn, V. Rivasseau and E. Trubowitz, An infinite volume 
expansion for many fermion green's functions. Helv. Phys. Acta65, 679-721 
(1992) 

[26] J. Magnen, G. Poirot and V. Rivasseau, The Anderson model as a matrix 
model, Nucl. Phys. B 58, 149-162, (1997) 

[27] Constructive Physics, Lecture Notes in Physics 446, Springer Verlag (1995) 

[28] V. Rivasseau, The two dimensional Hubbard Model at half-filling: I. Con- 
vergent Contributions, Journ. Stat. Phys. 106(3): 693-722; (2002) 

[29] G. Poirot, Mean Green's function of the Anderson model at weak disorder 
with an infrared cutoff. Ann. Inst. H. Poincare 70, 101-146, (1999) 

[30] J. Magncn, G. Poirot and V. Rivasseau Ward type identities for the 2d 
Anderson model at weak disorder, J. Stat. Phys. 93, 331-358, 1998 

[31] M. Disertori and V. Rivasseau, in preparation 



52 



